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Laminar Boundary-Layer Oscillations and 
Stability of Laminar Flow 


G. B. SCHUBAUER* ann H. K. SKRAMSTAD?t 
National Bureau of Standards 


ABSTRACT 


An account is given of an experimental investigation conducted 
at the National Bureau of Standards in 1940 and 1941, in which 
sinusoidal velocity fluctuations were discovered in the laminar 
boundary layer of a flat plate and in which the characteristics of 
these fluctuations were studied in detail and found to agree with 
the characteristics predicted earlier by the Tollmien-Schlichting 
stability theory. This work confirmed the theory, settling a con- 
troversy of many years and throwing new light on the causes of 
transition from laminar to turbulent flow. The work was done 
with the cooperation and financial assistance of the National 
Advisory Committee for Aeronautics and was published originally 
asan N.A.C.A. Advanced Confidential Report. The present re- 
view has been made possible by the declassification of the original 
report. 

The fluctuations are termed “laminar boundary-layer oscilla- 
tions” to distinguish them from the irregular velocity fluctuations 
previously observed by other investigators. In earlier investi- 
gations these oscillations were not found, probably because of the 
high turbulence in the wind tunnels in which investigations were 
conducted. The present investigation was conducted in a wind 
tunnel of turbulence much less than 0.1 per cent. The oscilla- 
tions were readily detected and studied by means of the hot-wire 
anemometer. 

A description is given of the methods used to produce and 
study boundary-layer oscillations. By these methods the oscil- 
lations are found to consist of a wave motion in the boundary 
layer. Amplified, damped, and neutral waves are found, and 
these Characteristics, together with wave velocities and type of 
Wave motion, are in accord with theory. Attention is called to 
the subsequent theoretical work of Lin, and some of Lin’s results 
are given. 


SyMBOLS 


* = distance from leading edge of flat plate 

y = distance from surface of flat plate 

Uy = mean velocity outside boundary layer 

4% = instantaneous x-component of fluctuation velocity 

0 = instantaneous y-component of fluctuation velocity 


Presented at the Aerodynamics Session, Fifteenth Annual 
Meeting, I.A.S., New York, January 28-30, 1947. Received 
August 30, 1946. 

— Aerodynamics Section, Division of Mechanics and 
nd. 

t Chief, Guided Missiles Section. 


= instantaneous fluctuation velocity perpendicular to Us 
and parallel to surface of flat plate 


root-mean-square values of u, v, and w 


= wave velocity 

= 2xf, where f = oscillation frequency 

= 27r/\, where \ = wave length 

= kinematic viscosity 

= boundary-layer thickness 

= boundary-layer displacement thickness 


8* = 0.3415 for Blasius velocity distribution 


= Uos*/v = boundary-layer Reynolds Number 
R. = Ugx/v = x-Reynolds Number 
R = 1.72+\/R; for Blasius velocity distribution 
q:/qo = ratio of dynamic pressure just outside boundary layer 
at any point x to dynamic pressure at arbitrary 
reference point. 


INTRODUCTION 


A pos STABILITY OF LAMINAR FLOW and the origin of 
turbulence has attracted widespread interest in 
recent years with the advent of the laminar-flow airfoil. 
However, long before laminar flow had assumed the 
practical importance it has today, transition from lami- 
nar to turbulent flow in boundary layers was of great 
scientific interest and was the subject of many investi- 
gations and a great deal of speculation. A number of 
investigators attempted to solve the theoretical prob- 
lem of the stability of laminar flow by determining un- 
der what conditions small disturbances in the form of 
velocity variations would increase or decrease with 
time. A brief account of such theories is given by 
Prandtl! with numerous references to works dating 
from that of Lord Rayleigh in 1880. This method of 
attack was recognized as the most logical one but was 
none too successful because of an oversimplification of 
the problem to avoid mathematical difficulties. 

About 17 years ago Tollmien? of Géttingen made im- 
portant advances in the application of the mathematical 
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theory and published a theory of the stability of laminar 
motion in the boundary layer near a thin flat plate in 
an air stream flowing parallel to the plate. According 
to this theory, small disturbances in velocity of any 
wave length lying within a certain region would be 
amplified, whereas disturbances of shorter or longer 
wave length would be damped. The calculations were 
repeated and extended by Schlichting* * in 1933 and 
1935. The amplified disturbances were assumed to 
grow until they caused a breakdown of the laminar flow. 
This was assumed to be the fundamental process of 
transition from the laminar to the turbulent state. The 
theory was universally discredited by scientists and 
engineers with a practical turn of mind, mainly because 
no evidence of disturbances with these characteristics 
had ever been found in boundary layers. All of the 
experimental evidence up to 1940 supported the view 
that any kind of disturbance, whether from turbulence 
in the air stream or irregularities in the surface, would 
cause transition when the disturbances were large 
enough or the Reynolds Number was sufficiently high. 

When investigators began to use hot-wire anemom- 
eters to study velocity fluctuations, they found irregular 
motions of large amplitude in laminar boundary layers. 
The reporting of this phenomenon by Dryden at the 


‘Fourth International Congress for Applied Mechanics 


in Cambridge, England, in 1934 attracted considerable 
attention. A discussion of’ the subject is given in 
reference 5. The fluctuations were recognized as dis- 
turbances impressed on the layer from the turbulence in 
the surrounding flow, and, while they eventually 
caused transition, they appeared to have no connection 
with the Tollmien theory. They were random and 
broke into turbulence with no apparent amplification. 
They were, in fact, additional evidence against the 
theory that transition resulted from the growth of small 
disturbances of a particular wave length, and they 
supported the growing conviction that stability or in- 
stability was merely a question of the size or intensity 
of the disturbance. 

The experimental investigation described here clari- 
fied this picture by the discovery of harmonic oscilla- 


tions in the boundary layer associated with waves wi nial 
properties that completely verified the theory of sele 
tive amplification and damping. The work was cop 
pleted late in 1941 and was published originally in refer 
ence 6. The present review has been made possible hy 
the declassification of the original report. The readg 
should consult reference 6 for a complete description 
the experiment and detailed results, as well as a revie 2 
of stability theory. After this experimental work ha 

been completed, Lin,’ studying at the California Insti 
tute of Technology, undertook a revision of the mathe 
matical theory of tlie stability of two-dimensional parz | 
lel flow and a clarification of some of the features of tly 
Tollmien-Schlichting theory which had been adverse} 
criticized. Some of Lin’s results are given here in Fig 


12. 
First EvIpENCE OF AMPLIFIED OSCILLATIONS 


Fic. 2. 


turbule: 
hot-wire 


In 1940 a research program was undertaken at th ramen 


National Bureau of Standards, with the cooperaticg; 


and financial assistance of the National Advisory Con . oot 
mittee for Aeronautics, to investigate the effectivenes irae 
1 


of damping screens in reducing wind-tunnel turbuleng 
and, at the same time, to investigate transition on a fla rere 
plate with different degrees of stream turbulence dom «<0 

to the lowest level attainable. The investigation wa Further 
conducted in the Bureau’s 4'/2-ft. wind tunnel show 
in Fig. 1. . A flat aluminum plate, !/, in. thick, 41/, 


screens, 


position 


wide, and 12 ft. long with a pointed leading edge (se " a 
Fig. 9), was placed vertically in the working chamberd cided to 
the tunnel. The pressure gradient along the directio in the le 
of mean flow was reduced to zero, on the average, }j sill on 
adjustable side walls near the tunnel walls. The posi neni 
use 
tion of transition was determined in the customary mai elite 
ner by a small pitot tube in contact with the surface d S 
ranged i1 
the plate by which the dynamic pressure was deter pe 
mined very near the surface at different distances fron ie A : 


the leading edge. The damping screens were placed i 
the settling chamber of the tunnel as indicated in Fig 
1, the number and type being varied over a wide range 
For each screen arrangement the turbulence of the fre 
stream opposite the plate was measured with the ust 
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Fic. 1. Elevation view of 4'/:-ft. wind tunnel. 
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Fic. 2. Effect of turbulence on x-Reynolds Number of transi- 
tion. Flat plate, zero pressure gradient. 


turbulence-measuring equipment, employing special 
hot-wire anemometers for obtaining all three compo- 
nents of the fluctuations. 

The results of these measurements are shown in Fig. 
2. As the turbulence of the free stream was reduced by 
increasing the number or the solidity of damping 
screens, the transition region moved progressively to 
larger values of the x-Reynolds Number until a value of 
the turbulence of about 0.08 per cent was reached. 
Further reduction in turbulence had no effect on the 
position of transition. 

When the stream turbulence had been reduced to a 
value in the neighborhood of 0.03 per cent, it was de- 
cided to observe what, if any, fluctuations were present 
in the laminar layer. It was felt that these would be 
much smaller than the fluctuations observed by Dryden 
because of the greatly reduced stream turbulence. For 
making this observation the hot-wire pickup was ar- 
ranged in the form of a “‘bug’’ as shown in Fig. 3. This 
arrangement slides fore and aft along the surface by 
remote control and places the hot wire at any desired 


Fic. 3. Hot-wire pickup for studying velocity fluctuations in 
boundary layer. 


value of x. A small celluloid block or sled, which rests 
directly on the surface, carries one or more wires at 
fixed distances from the surface on the forward extend- 
ing prongs. The wires were parallel to the surface and 
normal to the mean direction of flow so as to be sensi- 
tive to the u-component of the fluctuations. This 
arrangement of pickup was used to avoid relative mo- 
tion between the wire and the surface which might give 
false fluctuations due to vibration. The hot wires were 
connected to a compensated amplifier, just as for a 
turbulence measurement, and the output of the ampli- 
fier was connected to a cathode-ray oscillograph, as well 
as to the usual thermocouple and microammeter, so 
that the wave form of the fluctuations could be observed 
as well as the mean-square value. 

Irregular fluctuations were almost nonexistent as 
expected, but, as the pickup was moved downstream, 
an almost pure sinusoidal oscillation appeared, weak at 
first, but with increasing amplitude as the distance 
downstream increased. Just ahead of transition, 
bursts of extremely large amplitude occurred, and, at 
the initial point of transition, these bursts were accom-~- 
panied by a breaking into irregular high-frequency 
fluctuations characteristic of turbulence. These phe- 
nomena are illustrated in Fig. 4, which shows a set of 
film records made by photographing the oscillograph 
screen with a moving-film camera. Fig. 5 shows the 
oscillations at a lower wind speed and greater distance 
from the leading edge. Both figures show the progres- 
sion of events leading up to transition. 

Since the occurrence of the oscillations was unex- 
pected, a number of tests were made to make certain 
that the oscillations were not merely some effect of 
vibration. The possibility of such effects was soon 
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Fic. 4. Oscillograms showing laminar boundary-layer oscilla- 
tions in boundary layer of flat plate. Distance from surface = 
0.023 in.; Up = 80 ft. per sec.; time interval between dots = 
1/39 sec. 
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_ Fic. 5. Oscillograms showing laminar boundary-layer oscilla- 
tions in boundary layer of flat plate. Distance from surface = 
yo in.; Uy) = 53 ft. per sec.; time interval between dots = 

a sec. 
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Fic. 6. Effect of falling and rising pressure on laminar 
boundary-layer oscillations. Scale at upper left is ratio of pres- 
sure change to free-stream dynamic pressure. Distance from 
surface = 0.021 in.; Up = 95 ft. per sec.; time interval be- 
tween dots = 1/s9 sec. 


tuled out. In the course of the observations it was 
noted ,that the oscillations were definitely connected 
with transition, since the zone in which they occurred 
always preceded transition and moved with it fore and 
aft along the plate as the wind speed was varied. This 
was not absolute proof that the oscillations were the 
cause of transition, since there existed the possibility 
that the boundary layer became shock-excited by tran- 
sition occurring a short distance downstream, giving rise 
to an oscillation that was possibly the result of transi- 
tion rather than the cause. This was ruled out by re- 
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moving transition by means of an abrupt pressure drogimust f 
and yet leaving the oscillations totally unchanged at agiiposition 
upstream position. The results of this experiment anfifrequen 


shown in Fig. 6, which also serves to illustrate how @fharmon 
falling pressure damps out the oscillations. frequen 
All evidence pointed to the conclusion that thesiMparticul 


oscillations had some connection with the Tollmien{#pranch 


Schlichting stability theory. Accordingly, the theolfhave ut 
retical diagram shown in Fig. 7, taken from Schlichtfwhere | 
ing’s treatment,? was used as a test. This theoreticdamount: 
diagram applies to a Blasius velocity distribution ang Thus 
could be applied here, since the velocity distributiog§ count 
along the plate with zero average pressure gradient waaccordin 
found to be a Blasius distribution to within the accuragfnitial d 


of velocity measurements. As the oscillation frequenf,mount 
cies were determined from numerous oscillograms, sucifiictream. 
as those shown in Figs. 4 and 5, By /Uo? and the comiiicerved ir 
responding R were calculated and plotted on the diaMianswere 
gram. As shown in Fig. 7, the experimental points altyrbulen 
fell along Branch II of the curve. were fot 
The interpretation of this diagram is as folloy higher, t 
Within the loop ail disturbances are amplified, while iticayse of 1 
all other regions the disturbances are damped. Sincifappeared 
the curve divides the amplified from the damped, tl point tal 
curve itself is called the neutral curve. According tin good a 
theory, and later found to be true experimentally, 2 previous 
disturbances travel in a downstream direction with order of 
velocity given by their frequency and wave lengthiwas proc 
Small disturbances consisting of many frequencies (ranfMaminar 
dom disturbances, for example), originating in somselective 
upstream position of the boundary layer, all advance 1, the 
the direction of increasing R. As they advance, theffinterpret 
are first damped, then amplified as they traverse 
enclosed region, and then are again damped when they, ves we 
pass beyond Branch II. They are therefore most Ar 
highly amplified when they reach Branch II. A yiew to 
quency observation was always made at some partict§ pontrofte 
lar value of R. Taking the theory at face value, ClBang accic 
culations show that from a band of frequencies arriviny 
at a given value of R, the one most highly amplifie 


METHODS 


In the 
boundary 
completel 
using sou: 
the compl 
sound is 

wind tun 
The fir 
Speaker 
of variou 

the speak 


0 400 800 1200 1600 2000 2400 2800 3200 3600 backgrou 


Fic. 7. Theoretical neutral curve and amplification 2 lations ¢ 
according to Schlichting including experimental points {rom along thd 
frequency of oscillations on oscillograms. Theoretical to re 
amplification from Branch I to Branch II is indicated. quenc 
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LAMINAR BOUNDARY- 


ire dromtmust fall on Branch II. Therefore, because of the 
pd at almposition of the points, it is concluded that the observed 
hent angifrequencies were the ones most highly amplified. The 
p> how dfharmonic character of the records, or the purity of the 
frequency, may be accounted for by the fact that at a 
it thesfMparticular value of R all frequencies not falling on 
ollmienBranch II must lie either above the curve where they 
e theolfhave undergone some damping or below the curve 
chlichtRwhere they have not been fully amplified. This 
oreticdamounts to a filtering process. 
ion an Thus the observed oscillations could at last be fully 
ributiowligccounted for on theoretical grounds. They were, 
ent walaccording to theory, the amplified components of some 
initial disturbance, probably coming from the small 
amount of turbulence still remaining in the tunnel 
1S, Suciictream. The reason why such phenomena were ob- 
the coriiserved in this experiment and not in previous ones was 
she diafManswered by producing a tenfold increase in the stream 
urbulence and again looking for the oscillations. They 
were found, but, if the turbulence had been much 
followsiihigher, they would have been difficult to identify be- 
ause of the near coincidence of the point at which they 
SinofMappeared and the transition point. The one observed 
ed, thiEpoint taken at the higher turbulence is shown in Fig. 7 
ding tin good agreement with the other points. In all known 
ally, diprevious experiments the stream turbulence was of the 
order of 1 per cent. Under these conditions transition 
was produced by the initial disturbances, and the 
laminar layer was destroyed before any observable 
selective amplification could occur. 
In the experiment so far, the theory was used to 
interpret the results. It still remained to find out 
whether amplified neutral and damped oscillations or 
waves were a reality and whether the theory was cor- 
rect. An investigation was therefore begun with a 
view to setting up waves in the laminar layer under 
controlled conditions, rather than being left to random 
and accidental disturbances. 


METHODS OF PRODUCING OSCILLATIONS AND RESULTS 


In the search for schemes to excite oscillations in the 
boundary layer, a number of devices were tried before 
completely satisfactory results were obtained. Methods 
using sound, both pure notes and random noise, were 
none too satisfactory because of resonance effects and 
the complexity of the wave pattern in the tunnel. The 
experiments are worth mentioning, however, because 
sound is an unavoidable type of disturbance both in 
wind tunnels and in flight with power-driven aircraft. 

|_| The first scheme tried was to place a 25-watt loud- 
Speaker in the tunnel and “‘fill’” the tunnel with sounds 
of various intensity and frequency. The sound from 
|__JiBthe speaker was, of course, superimposed on the general 
noise of the tunnel. Boundary-layer oscil- 
x ote ations could be induced at will in a variety of positions 
's fo@along the plate by choosing the right combination of 
frequency and wind speed. Also, transition could be 
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DISTANCE BEHIND HOLE IN INCHES 
Fic. 8. Frequency of oscillations excited by sound through 
hole in plate. Up = 23 ft. per sec.; frequency expressed in cycles 
per second. 


moved 1 or 2 ft. ahead of its normal position. Random 
noise from the loud-speaker produced similar results, 
although the control was poorer. In general, the 
quantitative results were similar to those shown in Fig. 
7. 

The second scheme tried was an attempt to localize 
the sound field by bringing sound into the boundary 
layer through a small hole in the plate. A small head- 
phone was mounted on the back side of the plate oppo- 
site a '/s-in. hole. A survey of possible types of meas- 
urements showed that this method could be used only 
for obtaining the type of result given in Fig. 8. The 
hot-wire pickup was arranged to move upstream or 
downstream in the boundary layer downstream from 
the hole. The amplified output of the hot wire was fed 
back into the headphone, causing the system to oscil- 
late with a frequency depending on the distance between 
the hot wire and the hole, as shown in Fig. 8. Best re- 
sults were obtained for a wind speed of 40 ft. per sec. or 
less. The explanation of this performance is that 
proper phase relations had to exist for regeneration, and 
the frequency automatically decreased as the path 
length between the hot wire and hole increased in order 
to keep an integral number of waves between the 
source and the pickup. There existed a preferred fre- 
quency band, determined by the characteristics of the 
wire, electrical and acoustic circuits, and boundary 
layer. Jumps occurred when one additional wave 


length was necessary to keep the frequency near the , 


center of the preferred band. Reversing connections 
on the headphone changed the phase by 180° and re- 
sulted in the dashed curves. The distance between 
curves of any one type along any line of constant fre- 
quency is thus the wave length corresponding to that 
frequency. This showed for the first time the wave 
nature of the oscillations but otherwise conveyed little 
new information. 

Since the theory was known only for two-dimensional 
disturbances in a boundary-layer type of flow, what was 
really wanted was a source of two-dimensional disturb- 
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Fic. 9. Setup for exciting oscillations in laminar boundary 
layer by vibrating ribbon. 


ance. This was found in the vibrating-ribbon method 
illustrated in Fig. 9. A ribbon of phosphor bronze, 
0.002 in. thick and 0.1 in. wide, was mounted under 
tension on the plate as shown. Spacers were placed 
under the ribbon producing a 12-in. segment near the 
center uniformly parallel to the plate and about 0.006 
in. from the surface. When an alternating current was 
passed through the ribbon in the presence of a constant 
magnetic field from the electromagnets on the opposite 
side of the plate, the 12-in. segment vibrated to and 
from the surface in a single loop. The vibration set up 
a wave with the frequency of the ribbon which traveled 
downstream through the boundary layer and passed by 
the hot-wire pickup. The current through the ribbon 
was supplied from an oscillator, and the frequency 
could be varied through wide limits. The effect of the 
ribbon on the mean flow was exceedingly small, and 
what effects there were could not be detected 2 in. 
downstream. Neither the interferenee effects nor the 
results produced were critical to the dimensions and 
position of the ribbon. The ribbon was placed at 
various locations on the surface to cover the desired 
range of Reynolds Numbers. 

Two hot wires, one 0.010 in. and the other 0.110 in. 
from the surface, both sensitive only to the u-compo- 
nent, were mounted on the prongs of the “‘bug”’ so 


_ that it would be possible to study the wave proper- 
' ties at two fixed distances from the surface at vari- 
ous distances downstream from the ribbon along the 


centerline of the plate. Both hot wires were within the 


layer, and the velocity fluctuations introduced into the 
boundary layer could be picked up by either wire, 
Generally, only one wire was used at a time, depending 
on which gave the stronger response. Distance from 
the ribbon, wind speed, and frequency could be varied 
independently. The wave form showed some distor- 
tion for an inch or so behind the ribbon but became 
nearly sinusoidal at greater distances. Amplified, 
damped, and neutral waves were found at once, and it 
only remained to make quantitative measurements of 
amplitude, frequency, and wave length to obtain the 
necessary data for checking the theory. 

The usual procedure was to measure the root-meaa- 
square value of the u-component of the fluctuations, 
denoted by wu’, at various distances from the ribbon 
for some fixed wind speed Up at first one frequency, 
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Fic. 10. Growth and decay of u-component of oscillations 
produced by vibrating ribbon 4 ft. from leading edge of plate. 
xo is 2 in. downstream from ribbon; Uy = 64 ft. per sec.; fre 
quency expressed in cycles per second. 
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Fic. 11. Neutral curves (dashed) obtained from observed 
neutral frequencies compared to theoretical curve according to 
Schlichting (solid). 
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then another. In this way results similar to those illus- 
trated in Fig. 10 were obtained. In this figure, xo 
denotes a position 2 in. behind the ribbon and um’ is. the 
root-mean-square velocity fluctuation at x. The 
curve for 40 cycles per sec. shows damping throughout. 
As the frequency is increased, amplification takes place, 
reaches a maximum at 120 cycles per sec., and then de- 
creases until 180 cycles per sec. again shows damping. 


In more accurate terms, positive slopes at any value of 


x mean amplified waves, negative slopes mean damped 
waves, and zero slopes mean neutral waves. From 
sets of curves like these, not only could the neutral 
waves be found for outlining the amplification zone but 
amplification and damping coefficients could also be 
determined for the entire region. The reader is re- 
ferred to reference 6 for the complete results. The 
experimentally outlined neutral curves are shown in 
Fig. 11, again compared to Schlichting’s theoretical 
curve. 

The wave length and wave velocity were determined 
in the following simple manner: The input to the rib- 
bon from the oscillator was connected to one pair of 
plates of the cathode-ray oscillograph, and the output 
from the hot-wire pickup was connected to the opposite 
pair. A stationary Lissajou figure consisting of a 
single loop was obtained, since the frequencies of both 
input and output were the same. As the spacing be- 
tween the ribbon and the hot wire was changed, the 
phase change could easily be seen. Since a phase 
change of 180° could be detected accurately, the dis- 
tance intervals corresponding to one-half wave length 
were measured to determine the wave length. The 
wave velocity was then obtained by multiplying the 
wave length by the frequency. 

Wave lengths and wave velocities were measured for 
the neutral oscillations, and these furnished additional 
experimental data for comparison with theory. Fig. 12 
shows the neutral curves on a wave-length basis (aé* 
= 276*/\, where \ is the wave length). Recently, 
Lin’s theoretical curve became available from reference 
7, and this is shown in Fig. 12 in addition to Schlich- 
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Fic. 12. Neutral curves (dashed) obtained from observed 
neutral wave lengths compared to theoretical curves (solid). 
Closed circles—Branch I; open circles—Branch II. 
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Fic. 13. Neutral curves (dashed) obtained from observed 
velocity of neutral waves compared to theoretical curve of 
Schlichting (solid). 


ting’s curve. The experimental points are in a little 
better agreement with Lin’s curve than with Schlich- 
ting’s. On the whole the agreement between theory 
and experiment is good. 

Fig. 13 shows the neutral curves on a wave-velocity 
basis. Comparison is again made with Schlichting’s 
theoretical curve. It will be seen that the neutral 
waves are propagated downstream through the bound- 
ary layer, with a velocity definitely related to the 
free-stream velocity but only about one-third as 
fast. 

When these experiments were being performed, each 
check with theory was a stimulating experience. ‘There 
was nothing so unusual about setting up a wavy disturb- 
ance in the boundary layer, but finding that this wavi- 
ness really constituted a unique wave phenomenon 
with properties determined by the boundary-layer flow 
was out of the ordinary. Clearly, the wave was a 
motion superimposed on the mean flow, but what kind 
of motion was this? According to Schlichting’s theory 
in reference 4, when the longitudinal motions were in 
one direction near the surface they were in the opposite 
direction farther out. In other words, the boundary 
layer was not disturbed as a whole, but in any one sec- 
tion the oscillation consisted of speeds increasing in a 
certain interval of y while simultaneously decreasing in 
another. 

This theoretical prediction was accordingly put to the 
test by getting simultaneous records of the u-fluctua- 
tions from one hot wire at a fixed position near the sur- 
face and from another wire set at various distances from 
the surface in the same cross section of the layer. The 
records are shown in Fig. 14. The numbers on the left 
show the position of the movable wire, and the upper 
trace in each set shows the wave picked up by it. As 
the wire was moved outward, the wave amplitude de- 
creased to zero but remained in phase with the wave 
near the surface until it had disappeared. At 0.165 in. 
from the surface there is scarcely any perceptible oscil- 
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Fic. 14. Simultaneous records obtained with two hot wires 
located 1 ft. downstream from vibrating ribbon showing phase 
reversal in u-component of oscillations. Lower trace obtained 
from hot wire 0.055 in. from surface. Ribbon 3 ft. from leading 
edge; frequency 70 cycles per sec.; Up = 42 ft. per sec. 
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Fic. 15. Distribution of amplitude of oscillations across 
boundary layer. Solid curves are theoretical according to 
Schlichting; points and dashed curves are experimental 


lation. As the distance was further increased, the wave 


reappeared but now 180° out of phase with the wave 
near the surface. 

This was pleasingly in accordance with theory, not 
only qualitatively but quantitatively, as shown by Fig. 
15. Schlichting had calculated the distribution of 
amplitude for two positions on the neutral curve— 
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namely, R = 894 for Branch I and R = 2,070 for 
Branch II. The measurements shown in Fig. 15 
corresponded as nearly as possible to these two posi- 
tions. 

With the measurement of wave length, velocity, 
phase relationships, damping, and amplification, as 
shown by the u-component of the fluctuations, the 
physical characteristics of the waves were fairly well 
defined. The v-component, required by continuity, 
might also have been studied. However, it was not 
studied because of experimental difficulties and the 
belief that little additional information could have been 
gained by doing so. For the same reason, no attempt 
was made to observe the correlation between u and v 
predicted by theory. The general applicability of the 
theories of Tollmein, Schlichting, and Lin can no 
longer be doubted. For the latest and most compre- 
hensive treatment of the subject Lin’s work is recom- 
mended. 

An attempt was made to observe the transition into 
turbulent flow caused by oscillations produced by the 
vibrating ribbon, and a number of oscillograph records 
of this phenomenon are presented in reference 6. The 
results gave the general impression that transition 
occurred when the amplitude of the oscillations was 
sufficiently large. This goes back to the original idea 
held before 1940, with the important additional infor- 
mation—namely, that now we know that small out- 
side or initial disturbances are linked to the disturbance 
that finally causes transition through amplified boun- 
dary-layer oscillations. 


EFFECT OF PRESSURE GRADIENT ON BOUNDARY-LAYER 
OSCILLATIONS 


With but one exception, the foregoing results were all 
obtained with zero pressure gradient and a Blasius 
velocity distribution for which the theory was most 
complete. The one exception is shown in Fig. 6, where 
it is seen that a falling pressure decreases the oscilla- 
tions and an increasing pressure increases them. This 
effect was studied in more detail for oscillations pro- 
duced by the vibrating ribbon with pressures varying 
along the plate as shown in Fig. 16. Curve A shows 
the condition termed zero pressure gradient. Curves 
B and C’show moderate positive and negative gradients, 
respectively, and D and E show large gradients. The 
results on neutral oscillations obtained with pressure 
distributions A, B, and C are shown in Fig. 17. No 
theoretical curves are shown in this case, the curves 
labeled “‘zero pressure gradient” being the experimental 
curves taken from Fig. 11. Gradient C had no notice- 
able effect on the position of the neutral points, whereas 
gradient B expanded the region of amplification into 
the previously damped zone above Branch II. 

Corresponding figures for distributions D and E 
could not be obtained because of the large effect of 
these gradients on the oscillations. The general impres- 
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LAMINAR BOUNDARY-LAYER OSCILLATIONS 


sion given by observing the oscillations was that ampli- 
fication always occurred in gradient D regardless of the 
frequency and Reynolds Number, while damping al- 
ways occurred in gradient E. In other words, for the 
Reynolds Number range investigated with gradient D, 
the amplification zone occupied the whole region of a 
figure such as 17, and for gradient E, the amplification 
zone lay beyond the highest Reynolds Number reached 
(maximum Reynolds Number about 2,600). Suffi- 
ciently large pressure gradients therefore have a marked 
effect on the oscillations, negative gradients having a 
damping effect and positive or adverse gradients having 
an amplifying effect. This is in agreement with effects 
commonly observed on transition—namely, that a fall- 
ing pressure delays or prevents transition whereas a 
rising pressure brings about early transition. Thus we 
see, in part at least, why it is that laminar flow is main- 
tained over the region of falling pressure on an air- 
foil. 

Pressure gradient is only one of several factors that 
affect the stability of boundary layers. Curvature of 
the surface is another factor that was not touched upon 
but which has subsequently been investigated by Leip- 
mann? at the California Institute of Technology. Liep- 
mann also investigated the effect of small roughness 
elements on the surface and showed how these may ex- 
cite the boundary layer into oscillation. 


CONCLUDING REMARKS 


Perhaps a clearer physical picture of boundary-layer 
oscillations may be obtained by regarding a laminar 
layer as a medium in which a particular kind of wave 
motion is propagated, the propagation constants de- 
pending on the nature of the layer. This so-called 
medium has both positive and negative absorption co- 
efficients. The analogy to a medium should perhaps 
not be carried too far, but, if a boundary layer is com- 
pared to a water surface, it is easy to see how conditions 
would rarely be so “‘calm’”’ that even tiny ripples would 
not exist. Practically, then, some degree of wave 
motion will always be present. The important ques- 
tion is whether little ripples will grow to sizable 
waves, eventually producing transition to turbulent 
flow. 

There are, of course, those cases where the initial dis- 
turbances are so great—due, for example, to turbulence 
or surface roughness—that transition occurs at once. 
There are also those cases frequently met in practice 
where an adverse pressure gradient is sufficient to cause 
laminar separation followed almost immediately by a 
turbulent layer. Engineers are aware of these two 


cases and avoid them when attempting to preserve 
laminar flow. The technology of flight has therefore 
advanced to the stage where amplified waves play an 
important part in the preservation of laminar flow. It 
is now important to know where zones of amplification 
and damping lie and the magnitude of the amplification 
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and damping in these zones. As we have seen, the 
zones depend both on the boundary layer and the fre- 
quency or wave length of the disturbances. Both 


are quantities over which some control may be exer-: 


cised. 

So far the theory has been worked out only for the 
simpler cases. Furthermore, the present experimental 
work, which is limited to the flow along plane surfaces 
with zero, positive, and negative pressure gradients, is 
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little more than a good beginning toward an under- 
standing of the many flow problems involving the 
stability of boundary layers. 
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Noise Reduction in Aircraft 


H. WAYNE RUDMOSE? ano LEO L. BERANEKt 
Cruft Laboratory, Harvard University 


SUMMARY 


A complete treatment of the sources of noise in aircraft, the 
calculation or estimation of their intensities, and means for ap- 
plying acoustical treatment are given. The paper is divided into 
six principal parts as follows: I. Introduction; II. Sources of 
III. Variation of Noise Levels with Flight Conditions; 
IV. Sound Levels Inside Cabins; V. Procedure for Designing 
the Acoustical Installation; and VI. Miscellaneous Considera- 
tions. 

The principal sources of noise are the propellers, aerodynamic 
turbulence, the exhausts, engine vibration, and sources such as 
wind leaks, ventilators, loose ash trays, etc. The noise levels 
produced by the propellers in the lower frequency bands (35-150 
cycles per sec.) can be computed from flight data. Having the 
propeller noise, the shape of the remainder of the sound spectrum 


outside the airplane can be estimated. By taking into account . 


the acoustical properties of the bounding surfaces of the cabin 
undergoing treatment and the calculated noise spectrum out- 
side the airplane, the average sound level in the cabin is deter- 
mined. The method permits the detection of areas of the cabin 
which are weak acoustically and reveals ways to use the acousti- 
cal treatment efficiently. A practical example is carried through 
to illustrate the procedure. 


I. INTRODUCTION 


URING WoRLD War II, the Army Air Forces Air 
Technical Service Command and the Navy Bu- 
reau of Aeronautics sponsored at Harvard University, 
under funds provided by the Office of Scientific Re- 
search and Development, a program of research on 
methods for quieting airplanes.. This research led to 
the publication of a report entitled ‘‘Principles of Sound 
Control in Airplanes’! which was generally circulated 
throughout the industry. Since that date additional 
data have been obtained on more recent aircraft and 
the original concepts have been re-examined. A num- 
ber of the previous charts and procedures have been re- 
vised and are contained in this and in three other 
papers?’ * which, as a group, are intended to cover the 
various phases of the overall problem of airplane quiet- 
ing. 

In the second paper of reference 2 it was shown that 
an acceptable measure of the adequacy of acoustical 
treatment in airplanes is the percentage speech articu- 
lation, i.e., the number of words correctly understood 
by a passenger out of each one hundred words spoken 
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* This paper is based on work done for the Office of Scientific 
Research and Development under Contract OEMsr-658 with 
Harvard University. 

+ Now with Electro-Acoustic Consultants, Dallas, Tex. 

t Now Guggenheim Fellow jointly at Massachusetts Institute 
of Technology and Harvard University. 
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by another passenger. It was also shown that the per- 
centage speech articulation can be estimated if the 
sound levels in the three octave bands lying between 
600 and 4,800 cycles per sec. are known. A new figure 
of merit called the speech interference level obtained by 
averaging arithmetically the levels in these three bands 
was discussed. 

It is known® that the difference between desired 
speech interference levels and those actually measured 
in untreated airplanes is considerable. The design engi- 
neer is faced with resolving this difference at the time 
the original plans for an airplane are laid out, and he 
needs answers to the following questions: (a) How 
noisy will the airplane be with no sound treatment? 
(b) How much sound treatment must be added. to 
achieve the levels desired for proper speaking com- 
fort? (c) If the total weight of acoustical treatment is 
excessive, what possible alternatives in basic design are 
there for reducing the noise levels? To answer these 
questions one must consider (1) the sources of noise, 
(2) the manner in which the intensities and spectrum 
of these sources vary with changes in design and dif- 
ferent flight conditions, (3) the mechanism by which 
the noise is transmitted into a given cabin, and (4) 
how existing acoustical materials may be utilized to 
achieve maximum sound reduction for any specified 
added weight. These topics are discussed in this paper. 
Because of the great mass of data that has been ac- 
cumulated, only illustrative examples will be chosen to 
emphasize conclusions. In those instances where in- 
adequate data are available the tentativeness of con- 
clusions will be stated. 


II. Sources or NOISE 


Propeller Noise 

The propellers are the chief source of low-frequency 
noise on most airplanes. As a typical example of this 
fact, two analyses made in the pilot’s compartment of 
a commercial transport using a continuously variable 
sound frequency analyzer with a 5-cycle-wide band 
width are shown in Fig. 1. At the lower power condi- 
tion it is seen that the noise is principally that pro- 
duced by the first, second, third and fourth harmonics 
of the 60-cycle fundamental of the propeller-tip passage 
frequency. For greater powers the higher harmonics 
increase in intensity and the noise produced in the 
cabin becomes much louder. Two other typical ex- 
amples of the increase in overall level and distribution 
of levels as a function of frequency are seen in Figs. 2 
and 3 for two classes of military bombers. 
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Fic. 1. Analyses of noise levels in a transport using a 5-cycle wide band pass filter. Comparison of the lower with the upper graph 


shows the increase in propeller noise with increase in power and tip speed. 


Figs. 2 and 3 show that for a given number of decibels 
increase of sound levels in the lower two octave fre- 
quency bands (40-75 and 75-150 cycles per sec.) there 
is approximately the same increase in levels in the higher 
octave frequency bands. As the propeller-tip speeds 
approach 900-1,000 ft. per sec., the levels in the higher 
bands increase at a greater rate than they do in the two 
lower bands. Under normal cruising conditions pres- 
ent-day aircraft are operated with propeller-tip speeds 
not in excess of 800 ft. per sec. Hence, it appears that 
under the specific condition of normal cruising speeds 
it can be assumed that the entire spectrum shifts up or 
down by the same amount as does the sound level pro- 
duced by the propellers in the two lower octave fre- 
quency bands. 

Throughout this paper the portion of the noise spec- 
trum lying between 40 and 150 cycles per sec. and pro- 
duced principally by the propeller will be called ‘“‘Re- 
gion A.” 


Aerodynamic Noise 


Reference to Fig. 1 shows that at normal cruising 
speeds the spectrum of the noise in the higher frequency 
bands is of a continuous frequency type, i.e., all fre- 
quencies are represented. In these bands the noise is 


produced largely by vortexes generated by the propeller 
tips and by aerodynamic noises accompanying turbulent 
flow of the wind stream over the surfaces of the plane. 
In special cases some of the noise is contributed by the 
exhausts. For convenience, all of the noise pro- 
duced as a result of vortexes or turbulence is lumped 
together under the general title of aerodynamic 
noise. 
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Fic. 2. Octave band analyses of noise in a bomber as a function 
of flight conditions. 
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NOISE REDUCTION 


Octave band analyses of the noise produced under 
normal cruising conditions in a number of airplanes 
were taken and from the calculated attenuating proper- 
ties of the-cabin walls, the spectra of noise were com- 
puted outside the airplane. These spectra had almost 
the same shape (within a few decibels) as a function of 
frequency for all of the planes. This spectrum is 
shown in Fig. 4 as octave band levels relative to the 
levels in Region A. The data of Fig. 4 were confirmed 
by measuring the spectrum of noise around the outside 
of one airplane a few inches from the fuselage as the 
engine was revved up on the ground. The spectrum 
computed from data taken in flight and the spectrum 
measured on the ground were alike in the first four oc- 
tave bands. Above 600 cycles per sec., the levels meas- 
ured on the ground fell off with increasing frequency 
because of the absence of aerodynamic noise. 

For convenience in this paper, the eight octaves are 
divided into the several regions shown in Fig. 4. In 
Regions A and B the octave levels outside the plane are 
about equal; they increase rapidly in Region S; and 
fall off slightly in Region H. 


Exhaust Noise 


Extensive data show that aircraft noise produced by 
exhausts equipped with collector rings is negligible com- 
pared to that produced by the two types of noise dis- 
cussed above. This statement may not be true for 
ejector types of exhausts and evidence supporting this 
statement was collected on one military bomber equipped 
with ejector-type exhausts. Analyses were made of 
the noise in the flight-deck compartment of this bomber 
both with a 5-cycle-wide filter and with octave filters, 
and it was definitely found that the exhaust noise ex- 
ceeded the propeller noise for cases when the propeller- 
tip speed was less than 700 ft. per sec. and the engine 
horsepower was less than half of take-off power. The 
data from these measurements are shown in Fig. 5. 

In Fig. 5 the measured noise in the airplane in the 
lower octaves is shown compared to the computed noise 
for the propellers and exhaust alone. At the lowest 
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tion of flight conditions. Note that the spectra are essentially 
parallel at the lower horsepowers. 
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Fic. 5. Octave sound levels in Region A vs. engine power show- 
ing relative intensities of propeller and ejector exhaust noise. 


horsepower, the measured levels are 6 db. higher than 
those computed from information given later in this 
paper as being produced by the propeller alone. The 
association of the levels in these low-frequency bands 
with components of exhaust noise was made possible 
by the 5-cycle-wide band frequency analyses. At about 
600 hp. the exhaust and propeller noise components be- 
came equal in level and above that point the propeller 
components were found to predominate. 

It appears likely from this limited evidence that the 
use of ejector exhausts may result in some cases in in- 
creased noise levels at normal cruising speeds. Also, 
because the exhausts are pointing backward, it is ex- 
pected that their influence will be greatest in the rear 
compartments of the plane. Data were not taken in 
the rear compartments on this particular airplane. 


Vibration Noise 


Inadequate engine mounts will result in the produc- 
tion of noise in the airplane cabin because of vibration 
of the fuselage. In at least one instance studied, an 
A-20 attack bomber, this transmitted vibration pro- 
duced noise levels in excess of those produced by the 
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propeller. The manner in which such noise can be 
transmitted into the cabin will be discussed in Section 
IV. 


Secondary Sources of Noise 


Miscellaneous sources of noise include rattles pro- 
duced by vibrating ash trays, magazine racks, etc.; 
wind leaks around doors, windows, and turrets; noise 
created by the ventilating system; and noise created 
by electrical generators, air compressors, etc. In gen- 
eral, each of these sources of noise must be quieted in- 
dividually.to such a degree that its contribution to the 
total noise is less than that of the noise entering the air- 
plane through the acoustically treated sidewalls. Fail- 
ure to pay due attention to such small details often re- 
sults in a nullification of a large part of the effectiveness 
of an acoustical lining, and generally the cost in terms 
of weight needed to quiet these secondary sources is 
negligibly small in comparison with the weight of treat- 
ment required to reduce the noise originating from the 


propellers. 
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Fic. 7. Effect of eliminating wind leaks around windows on the 
noise inside the cabin of a XPB2Y-3. 
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Examples of instances where failure to attend to de. 
tails resulted in substantially increased noise levels are 
shown in Figs. 6 through 9. 

In the example of Fig. 6, the rear doorway was poorly 
fitted and the wind generated a loud hissing noise. In 
the washroom on the opposite side of the plane, the 
noise levels in the higher octave bands were as much as 
18 db. lower than those produced near the door. 
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Fic. 8. Comparison of noise levels near and remote from the 
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Fic. 9. Comparison of noise levels at the pilot’s position in a 
PBM.-3 with ventilator open and closed. 


Window leaks can cause similar results as is demon- 
strated in Fig. 7. In this plane two escape hatches were 
above the heads of the pilots. These hatches opened 
outward and were held down by ordinary twist locks. 
During flight, the reduced pressure on the top side of 
the airplane raised these hatches far enough to open a 
1/,-in. air leak around the edges. The high-frequency 
hissing noise that resulted raised the noise levels by 
almost 25 db. in the higher frequency octaves. 

Wind leaks around turrets produce similar effects, 
as can be seen for the case of Fig. 8. In this instance, 
the sound levels were measured in the same compart- 
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ment at two positions, one near the turret and the other 
somewhat removed. 

Ventilators, unless properly designed, can be exces- 
sively noisy. Data are shown in Fig. 9* for a bomber, 
in which case the noise produced by the ventilator ex- 
ceeded by several decibels that produced in the higher 
frequency octaves by all other sources of noise. 


VARIATION OF NoIsE LEVELS WITH FLIGHT 
CONDITIONS 


III. 


Sound Levels vs. Engine Power and Propeller-Tip Speed 


It has already been pointed out that, for normal 
cruising conditions, the levels outside the airplane in 
Regions B, S, and H are related to those in Region A 
according to the curve shown in Fig. 4, and that the 
levels in Region A are determined in most cases by the 
propeller noise. Hence, it is important that relations 


10 T T 


Sound Level vs. Vo 
for a Constant 
Power of 600 HP. 


XBP2yY-3 


™ 600 700 800 900 
PROPELLER TIP. SPD. 
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Fic. 10. Contours of noise levels in Region A vs. propeller tip 
speed with engine power held constant at 600 hp. 


between propeller noise, flight conditions, and location 
of the propeller with respect to the cabin under con- 
sideration be known. 

During the course of the war, measurements of noise 
levels were made in about 50 different types of airplanes 
under a wide variety of flight conditions to obtain such 
telations. Sufficient data were obtained in only seven 
types to permit a detailed determination of the pa- 
rameters affecting noise levels in Region A. 

In correlating the datum points for these different 
planes with each other, it is first necessary to set down 
a list of conditions under which comparisons of data 
may be made. These conditions are: 


* The authors are indebted to H. L. Ericson and R. L. Wallace, 
Jr., for the data shown in this figure. 
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_ Fic. 11, Contours of noise levels in Region A vs. engine power 
with propeller tip speed held constant at 500 ft. per sec. Note 
the different propeller tip separations. 


(a) Three-blade propellers. 

(b) Two or four engines. 

(c) Collector ring exhausts. 

(d)_ Microphone positions within 6 ft. of the plane 
of the propellers and about 2 ft. from the fuselage wall 
at the head level of a sitting man. 

(e) Reasonable seals around doors, windows, and 
hatches. 

(f) Altitudes less than 10,000 ft. 

(g) Good vibration insulation of the engines from 
the wings. 

Modification of the conclusions found under the re- 
strictions mentioned above were studied as separate 
items and are treated later in this paper. 

Contours of sound levels vs. power developed per 
engine with propeller-tip speed held constant were ob- 
tained on three of the planes at a position in each plane 
corresponding to (d) above. From these contours the 
curves of Fig. 10 were plotted for the condition of a con- 
stant power of 600 hp. with the propeller-tip speed as 
the independent variable. These contours lie parallel 
to each other, showing that the sound level in Region 
A (i.e., the arithmetic average of the levels in the first 
two octave bands) increases approximately 2.7 db. for 
each 100 ft. per sec. increase in propeller-tip speed. 
This fact was then used to correct the data on seven 
planes to an equivalent propeller-tip speed of 500 ft. 
per sec. by subtracting from each datum point a num- 
ber of decibels equal to 2.7(Vo — 500)/100 where Vo 


Vo* 500 FT. PER SEC. 

= 

8 INCHES 


AVERAGE OCTAVE LEVEL IN 08 


POWER PER ENGINE IN HORSEPOWER 


Fic. 12. Contours of noise levels in Region A vs. engine power 
with propeller tip speed held constant at 500 ft. per sec. Data 
are corrected to a tip separation of 8 in. 
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Fic. 13. Computational charts for determining noise levels in Region A vs. engine power, propeller tip speed, and minimum pro- 
peller tip clearance from the fuselage. 
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is the tip speed for that point. These results are plotted 
in Fig. 11. For all the airplanes, the sound levels in- 
creased at a rate of about 5.5 db. for each doubling of 
horsepower. The curves do not lie on top of each other, 
but if the minimum clearance of the propeller tips from 
the fuselage is taken to be a major factor in influencing 
the noise levels, the curves can almost be made to coin- 
cide. 


Sound Levels vs. Minimum Tip Clearance 


The curve relating, d, the minimum propeller-tip 
clearance, to the correction in sound level needed to 
bring the curves into coincidence is shown in Fig. 13b. 
Many datum points additional to those for the airplanes 
of Fig. 11 were selected from data taken on the 15 or so 
additionally comparable airplanes studied and used to 
determine the exact shape of this curve. The scatter 
of the points on this graph was large and is probably 
the result of differences in propeller-tip shapes, the ef- 
fect of standing waves in cabins, etc. The data of Fig. 
11 are shown corrected in Fig. 12 according to the curve 
of Fig. 13b and are in remarkable agreement. 

The conclusions reached from this analysis are that 
the noise produced in an airplane cabin in Region A is 
generated primarily by the propeller tip that is nearest 
the fuselage at any given instant. Otherwise, if this 
were not true, the levels would not decrease rapidly as 
the minimum propeller-tip separation is increased. 
When the separation becomes large enough, as would be 
expected, the noise appears to come from a larger and 
larger source, i.e., more of the arc traced by the pro- 
peller tips, and the levels drop off much less rapidly. 
From Fig. 13b we can conclude that the desired mini- 
mum separation of the propeller tips from the fuselage 
is of the order of 20 in.; and that separations greater 
than this are not necessary from an acoustical stand- 
point. 

The complete formula for calculating the sound levels 
in Region A under the assumptions previously made is 
as follows: 


43 24 HP 
S.L. = 102 + y or + 18.3 log 800 + 
2.7 (Vo — 700) 

where 

S.L. = average octave sound level in Region A 

dy = minimum propeller-tip clearance 

HP = power delivered by each engine in horsepower 

Vo = propeller-tip speed in ft. per sec. 


This equation in graphical form is presented in Fig. 13 
and forms the basis of a method for calculating sound 
levels in airplanes in advance of construction. 


Effect of Propeller-Tip Shape on Sound Levels 


If the statement of the previous paragraph is true, 
namely, that the sound levels in Region A are generated 


a 
= 
zo 

° 

- 

* 

c 

oO 


“80 -6€0 -40 -20 20 40 80 100 


VARIATION OF AIR SPEED IN M.P. H. 
FROM NORMAL CRUISING 


Fic. 14. Variation of noise levels in Regions S and H vs. varia- 
tion of air speed in m.p.h. from normal cruising. 


primarily by the propeller tip nearest the fuselage, then 
it is to be expected that the shape of the propeller tip 
will have some effect on the noise levels. The meager 
data at our disposal show that a wide propeller blade 
having broad tips produces more noise, by several deci- 
bels, than a narrow blade having finer tapered tips. 
Further study on this important point is indicated. 


Three- vs. Two- or Four-Blade Propellers 


If the conclusion is correct, that it is the noise pro- 
duced by the propeller tip nearest the fuselage which 
is important, then the charts of Fig. 13 ought to be 
valid for propellers with different numbers of blades 
provided that the numbers on the abscissa are multi- 
plied by a fraction equal to the ratio of three to the 
number of blades on the propeller under consideration. 
Hence, for a given engine power and tip speed, the 
levels for a four-blade propeller should be 2 db. less 
and for a two-blade propeller 3 db. more than those 
shown in Fig. 13a. Insufficient data were taken to de- 
termine how accurate this statement is but the data 
that have been taken indicate that these values are 
approximately correct. 
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Fic. 15. Variation in sound levels at 35,000 ft. from values 
obtained at 5,000 ft. Shaded region shows spread of data ob- 
tained at all crew positions in B-17F 
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is the tip speed for that point. These results are plotted 
in Fig. 11.. For all the airplanes, the sound levels in- 
creased at a rate of about 5.5 db. for each doubling of 
horsepower. The curves do not lie on top of each other, 
but if the minimum clearance of the propeller tips from 
the fuselage is taken to be a major factor in influencing 
the noise levels, the curves can almost be made to coin- 
cide. 


Sound Levels vs. Minimum Tip Clearance 


The curve relating, d, the minimum propeller-tip 
clearance, to the correction in sound level needed to 
bring the curves into coincidence is shown in Fig. 13b. 
Many datum points additional to those for the airplanes 
of Fig. 11 were selected from data taken on the 15 or so 
additionally comparable airplanes studied and used to 
determine the exact shape of this curve. The scatter 
of the points on this graph was large and is probably 
the result of differences in propeller-tip shapes, the ef- 
fect of standing waves in cabins, etc. The data of Fig. 
11 are shown corrected in Fig. 12 according to the curve 
of Fig. 13b and are in remarkable agreement. 

The conclusions reached from this analysis are that 
the noise produced in an airplane cabin in Region A is 
generated primarily by the propeller tip that is nearest 
the fuselage at any given instant. Otherwise, if this 
were not true, the levels would not decrease rapidly as 
the minimum propeller-tip separation is increased. 
When the separation becomes large enough, as would be 
expected, the noise appears to come from a larger and 
larger source, i.e., more of the arc traced by the pro- 
peller tips, and the levels drop off much less rapidly. 
From Fig. 13b we can conclude that the desired mini- 
mum separation of the propeller tips from the fuselage 
is of the order of 20 in.; and that separations greater 
than this are not necessary from an acoustical stand- 
point. 

The complete formula for calculating the sound levels 
in Region A under the assumptions previously made is 
as follows: 


HP 
S.L. 24 18.3 log 600 


de 
2.7 (Vo — 700) 
where 
S.L. = average octave sound level in Region A 
dj = minimum propeller-tip clearance 
HP = power delivered by each engine in horsepower 
Vo = propeller-tip speed in ft. per sec. 


This equation in graphical form is presented in Fig. 13 
and forms the basis of a method for calculating sound 
levels in airplanes in advance of construction. 


Effect of Propeller-Tip Shape on Sound Levels 


If the statement of the previous paragraph is true, 
namely, that the sound levels in Region A are generated 
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Fic. 14. Variation of noise levels in Regions S and H vs, varia- 
tion of air speed in m.p.h. from normal cruising. 


primarily by the propeller tip nearest the fuselage, then 
it is to be expected that the shape of the propeller tip 
will have some effect on the noise levels. The meager 
data at our disposal show that a wide propeller blade 
having broad tips produces more noise, by several deci- 
bels, than a narrow blade having finer tapered tips. 
Further study on this important point is indicated. 


Three- vs. Two- or Four-Blade Propellers 


If the conclusion is correct, that it is the noise pro- 
duced by the propeller tip nearest the fuselage which 
is important, then the charts of Fig. 13 ought to be 
valid for propellers with different numbers of blades 
provided that the numbers on the abscissa are multi- 
plied by a fraction equal to the ratio of three to the 
number of blades on the propeller under consideration. 
Hence, for a given engine power and tip speed, the 
levels for a four-blade propeller should be 2 db. less 
and for a two-blade propeller 3 db. more than those 
shown in Fig. 13a. Insufficient data were taken to de- 
termine how accurate this statement is but the data 
that have been taken indicate that these values are 
approximately correct. 
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Variation of Sound Levels-with Air Speed 


Data on about six airplanes have been taken on the 
variation of sound levels with indicated air speed. To 
obtain these data, a level flight condition for the air- 
plane was first established and the manifold pressure, 
engine r.p.m. indicated air speed, and the octave sound 
levels were recorded. Then the airplane was placed in 
a climb, maintaining the same manifold pressure and 
engine r.p.m.’s, but with reduced air speed. The plane 
was made to climb at different rates, and for each the 
indicated air speed and the octave sound levels were re- 
corded. Finally, the plane was placed in a dive, keep- 
ing the manifold pressure and engine r.p.m. constant, 
and the octave sound levels as a function of indicated 
air speed were recorded. 

The cumulative data show that the sound levels in 
Region A manifest a tendency to shift upward with 
indicated air speed, while for Regions S and H there is 
a definite shift of about 6 db. upward for each 100- 
m.p.h. increase. The data are summarized in Fig. 14* 
as corrections in decibels vs. indicated air speed in 
m.p.h. relative to the normal cruising speed for that 
particular airplane. To use this chart in estimating 
sound levels in airplanes, calculate the level in Region 
A from Fig. 13. Then correct the value obtained ac- 
cording to Fig. 14 if the plane is in a diving or climbing 
attitude. No other correction for air speed has been 
found necessary. 


Variation of Sound Levels with Altitude 
Extensive data were taken in a Boeing Fortress at 


nearly all of the crew positions for altitudes of 5,000 and ~ 


35,000 ft. A summary of the results is shown in Fig. 
15.* The graph shows the spread of results obtained 
by plotting the octave sound levels at 35,000 ft. relative 
to the levels at'5,000 ft. in each of the frequency bands. 
The plane was operated at 63 per cent rated power at 
5,000 ft. and at approximately 75 per cent rated power 
at 35,000 ft. For comparison at equal powers at these 
two altitudes, the envelope in Fig. 15 should be shifted 
downward by approximately 3 db. These results indi- 
cate that for the same horsepower delivered to the pro- 
peller, the noise is little different at 35,000 ft. than it is 
at low altitudes. 


IV. Sounp LEvEts INSIDE CABINS 


Theoretical Considerations: Air-Borne Sound 


In the previous sections evidence was given that the 
shape of the spectrum outside the airplane is about the 
same for all airplanes and that the level in Region A can 
be calculated by the use of Fig. 13. It now remains to 
determine the effect of the sidewalls and the windows 
of the cabin in reducing the noise. 


* The authors are indebted to H. F. Dienel for the data of Figs. 
14 and 15. 


It is at present not feasible to calculate sound levels 
in a cabin in which whistling noises are generated by air 
leaks to the outside. Also, it is difficult to calculate the 
levels in a cabin in which some areas are much less ef- 
fective in reducing transmitted sound than are others. 
Under the special condition, however, that the cabin is 
tightly sealed and that all areas of the bounding sur- 
faces of the cabin are nearly equally effective in reduc- 
ing transmitted sound, quite satisfactory calculations 
of the levels of sound in the cabin can be made. This 
condition of uniformity of sidewalls in terms of their 
acoustic properties is also the most efficient design. 
Psychologically, the ear is aware only of the loudest 
sound. impinging on it, and if the listener is near an 
acoustically weak area little good is to be derived from 
added treatment over other portions of the surface of 
acompartment. Hence, under this condition, inefficient 
use has been made of the added material. 

From fundamental acoustics a suitable formula is 
available for computing the difference between the 
sound levels inside and outside an acoustically uniform 
enclosure when this enclosure is immersed in a uniform 
sound field. The relation is: 


Db. sound reduction = 10 logy [1 + (a/7)] (1) 


where a and 7 are the average sound absorption and 
transmission coefficients, respectively, of the bounding 
surfaces of the cabin and are determined by the formu- 
las: 
a = (Sija’ + Spa” + + ....)/S (2) 
= (Sir! + Sor” + + ....)/S (3) 


In these formulas S,, S:, S;, etc., are the areas of por- 
tions of the boundaries having the same construction, 
and a’, a", r’, 7”, etc., are the particular a’s and 7r’s 
for these areas, respectively. Also, 


S = S, + S:+ S3 + .... = total boundary area (4) 


A graph of Eq. (1) is given in Fig. 16 for various values 
of rand a. If we let 


A = Saand T = Sr (5) 
then Eq. (1) becomes 
Db. sound reduction = 10 logy [1 — (A/T)] (6) 


where A is the total number of absorption units in the 
cabin and T is the total number of transmission units 
of the bounding surfaces. 

One other important restriction also applies to Eq. 
(6), namely, the basic derivation assumes that the 
acoustic energy arriving inside the cabin is diffused uni- 
formly throughout the volume. In actual airplane 
cabins, however, standing waves exist with a resultant 
nonuniformity of sound distribution? (see Fig. 20). 
These restrictions are important in that they serve a 
a warning that computations made using Eq. (6) must 
be interpreted properly. 


Fic. 


4 
7, 
| 


levels 
Dy air 
te the 
ss ef- 
thers. 
bin is 
sur- 
educ- 
itions 

This 
their 
esign. 
udest 
ir an 
from 
ice of 
icient 


ila is 
1 the 
iform 
iform 


NOISE REDUCTION IN AIRCRAFT 


ATTENUATION OF CABIN WALLS IN DB 


A 
VA 
A 
7 il 
| 26 
“e 28 
34 
38 
44 
= 
4 


ABSORPTION COEFFICIENT. Oc 


Fic. 16. Theoretical reduction of sound levels vs. absorption coefficient and either transmission coefficient (right-hand ordinate), 
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Agreement Between Measurement and Theory 


A detailed tabulation was made of the properties of 
the bounding surfaces in each of about 20 airplanes in 
which sound levels had been obtained. The ratio of A 
to T for several frequencies throughout the range was 


then computed. Using Eq. (6), the sound spectrum in- _ 


side the airplane was computed assuming the spectrum 
of Fig. 4 as existing outside the airplane. (The spectrum 
of Fig. 4 was determined by the inverse of this process 
for several of these 20 airplanes.) The agreement was 
satisfactory within reasonable tolerances (about +3 
db. in Region A, +5 db. in Region S, and +6 db. in 
Region H) provided the computations for each octave 
band were made using 7’s and a’s measured at the lower 
frequency limit of the octave band. For example, to 
compute the measured sound level in the 300-600 
cycles-per-sec. octave band, the values of r and a neces- 
sary were those measured at 300 cycles per sec. 

If acoustical treatments with surface densities less 
than 0.5 Ib. per sq.ft. are employed, the sound level in 
Region A remains relatively unaffected by the addition 
of material. However, if the transmission coefficient 
7 becomes appreciably less than unity and the absorp- 
tion coefficient a increases above zero at frequencies 
below 150 cycles per sec., then some reduction in level 
of Region A may be expected according to the chart of 
Fig. 16. 

Comparisons between calculated and measured sound 
levels in the eight octave bands between 40 and 10,000 
cycles per sec. for two transports and two bombers are 
shown in Fig. 17. Similar agreement between calcu- 
lated and measured levels was obtained for most other 
airplanes provided ventilator noises and noises caused 
by air leaks around hatches and turrets were taken into 
account. A tabulation of estimated vs. measured levels 
in Region A for twelve additional airplanes is made in 
Table 1. 


Theoretical Considerations: Structure-Borne Sound 


In the preceding sections, the tacit assumption has 
been made that the vibration that excites the skin is 
air-borne vibration. The values of 7 are taken from ex- 
perimental data obtained by applying a sound pressure 
to the primary side of the panel, the bounding edges 
being rigidly clamped. In an airplane, the panel is 
excited not only by air-borne sound but also by 
vibration that reaches the panel after traveling 
along the skin of the airplane, i.e., structure-borne vi- 
bration. 

To demonstrate the fact, assume that the principle 
source of the structure-borne vibration is the engine. 
Such vibration can be due to sharp impacts that occur 
in the operation of cams, tappets, etc. Also, engine un- 
balance, flexural and torsional vibrations, etc., are addi- 
tional sources. Let it be assumed that the spectrum of 
the source of structure-borne vibration has uniform 
energy per cycle (increases 3 db. per octave). Any 
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Fic. 17. Calculated sound spectra compared with sound. spec- 
tra measured in flight in the pilot’s compartment of four airplanes. 


sharp click or the sum of a number of sharp clicks would 
satisfy this energy requirement. 

In traveling from the source to the pertinent panel, 
these vibrations are attenuated by the elastic mounts 
used to support the engine. For effective isolation of 
the engine vibration from the remaining structure of 
the airplane, the resonance frequency of the engine and 
its suspension is chosen low, i.e., less than 50 cycles per 
sec. (3,000 r.p.m.). Hence, for vibrations having fre- 
quencies above this resonant frequency, the attenua- 
tion increases uniformly with increasing frequency (6 db. 
per octave increase). Such an attenuation corresponds 
to that suffered by air-borne vibration in traveling 


TABLE 1 
Meas- Esti- Differ- 

Vo, mated ence 

H Ft. do, Level Level in 

Plane Engine _ per Sec. In b. Db. 
PBM-1 750 725 17.5 104 105 +1 
DC-3B 600 740 8 106 108 +2 
980 118 115 -3 

900 810 24 107 107 0 

XPBY-5A 800 800 10 111 113 +2 
B-18 . 636 770 8.5 109 109 0 
B-17C 600 10 102 102 0 
640 595 10 104 104 0 

680 634 10 105 105 0 

705 6 10 106 106 0 

730 10 106 107 +1 

760 7 10 108 107 -1 

770 718 10 106 108 +2 

780 7 10 107 108 +1 

790 736 10 108 109 +1 

B-314 750 825 24 107 107 0 
C-53 695 828 8 112 111 -i 
828 8 110 110 0 

720 812 8 112 111 -1 

595 812 8 110 110 0 

725 778 8 110 110 0 

580 778 8 109 109 0 

705 745 8 109 109 0 

560 745 8 106 108 +2 

535 710 8 103 106 +3 

675 710 8 104 108 +4 

405 676 8 101 103 +2 

595 676 8 104 106 +2 

C-47 8 608 8 100 103 +3 
592 642 8 104 105 +1 

568 676 8 104 1 +2 

659 710 8 106 107 +1 

679 745 8 104 109 +5 

B-34 650 497 4.5 107 105 -2 
650 497 4.5 108 105 -3 

1,163 747 4.5 114 115 +1 

PBM-3 750 640 14 104 104 0 
965 675 14 106 107 1 

1,500 4 108 114 6 

XPBY-1 425 682 55/s 111 106 -5 
625 865 55/s 119 113 -6 

890 55/s 122 120 -2 

Vega 715 520 1.75 114 110 -4 
PVv-1 495 4.75 111 105 -6 
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through a dural panel when the frequencies of vibration 
are well above the first two or three resonant frequencies 
of the panel. This correlation is to be used later in the 
discussion. 

Continuing to follow the structure-borne path, the 
vibrations, after leaving the engine mounts, travel 
along the wing surface and then along the skin until 
they reach the relevant panel. Assume the worst case 
first, namely, that over this transmission path there is 
essentially no attentuation. Having reached the panel 
under consideration, the structure-borne vibration is 
radiated in the form of sound waves. The radiation 
traveling into the cabin suffers no attenuation as a 
result of the air cavity between the dural and the acous- 
tical material (assuming the acoustical material is 
spaced out from the dural). Thus the structure-borne 
vibration, once it leaves the dural in the form of radi- 
ated sound energy is attenuated only by the acoustical 
material before it reaches the interior of the cabin. 
Air-borne sound would be attenuated not only by the 
acoustical material but also by the dural panel; hence, 
the structure-borne vibrations have an easy path into 
the cabin by comparison with the path of air-borne vi- 
brations. 


On the basis of this latter argument, it might appear 
that structure-borne vibration should play a prominent 
role in determining the amount of sound energy which 
appears inside the airplane. However, no assignment of 
relative amplitudes of the two vibrations has been 
made and, furthermore, the structure-borne vibration 
is attenuated by the engine mounts before the vibration 
reaches the panel and in the actual case is attenuated 
in being transmitted from the mounts to the panel in 
question. Assume a severe case of vibration—one 
where for a frequency of 50 cycles per sec. the ampli- 
tude of the structure-borne vibration as it leaves the 
inside of the dural panel is equal to the amplitude of 
the air-borne vibration which strikes the outside of the 
dural. Such an assumption would mean that the spec- 
tra, expressed in terms-of energy per cycle, inside the 
cabin due to the two types of vibration would be that 
shown in Fig. 18. The curves cross at 100 cycles per 
sec. so that the structure-borne levels become increas- 
ingly greater below that frequency. 


The results shown in Fig. 18 naturally suffer in detail 
because they are “‘smoothed”’ curves, that is, panel 
resonances are averaged; however, the principle con- 
clusions are reasonable. It is noted that structure- 
borne vibration would tend to increase the sound levels 
in the two octaves below 150 cycles per sec. In fact, 
the lowest octave would be increased more than the 
75-150 cycles-per-sec. octave, making the levels in 
these two octaves unequal; hence, if the “below 75” 
is considerably higher than the 75-150 cycles-per-sec. 
octave it is probably reasonable to conclude that the 
structure-borne vibration cannot be neglected. The ef- 
fect of structure-borne vibration would be to raise the 
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Fic. 18. Theoretical calculation of air-borne and structure- 
Lage sounds arriving inside cabin under assumptions outlined 
in the text. 


spectrum above what would be expected on the basis of 
air-borne vibration. Thus, the calculated spectrum, 
based on the charts presented in this paper, would lie 
below the measured spectrum; the greatest discrepancy 
existing for the low-frequency octaves. It is believed 
that the spectrum of the A-20-B obtained by this labora- 
tory in 1941 is of this type although detailed vibration 
measurements were not taken at the time to confirm 
the conclusion (see Fig. 19). The structure-borne 
spectrum will drop off more rapidly than indicated in 
Fig. 18 for frequencies above 1,000 cycles per sec. be- 
cause some types of dissipation have been neglected in 
the structure-borne spectrum, whereas all measurable 
dissipation is included in the air-borne spectrum. 

The role that structure-borne vibration plays in pro- 
ducing noise inside the airplane is qualitatively ex- 
plained by this discussion. However, at present no 
scheme exists for estimating quantitatively the amount 
of structure-borne vibration which will exist in aircraft. 
If such a procedure is developed, the knowledge can be 
incorporated in the overall method of estimating sound 


spectra. 


130 oA. T 
3 A-20-B 
HAVOC 7 
z 100 
Ss 
a3 
~ 
CALCULA TED | |_ 
2 
2 
wt 
60 
50 
3 
150 300 600 1200 2400 f>4800 
4800 


75 
150 300 600 1200 2400 
FREQUENCY OCTAVES IN C.PS. 


Fic. 19. Measured levels in A-20-B attack bomber vs. cal- 
culated levels assuming no vibration. The difference between 


measurement and calculation is attributed to noise produced by 
engine vibration. 
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Fic. 20. Contours of sound levels in two octave bands in a DC- 
B transport. 


Sound Levels oi Position 


The sound levels in a compartment vary somewhat 
as a function of position. This is because of both the 
presence of standing waves and the location of the plane 
of the propellers. Typical data are shown in Fig. 20 
for sound levels in a DC-3B in regular air-line service. 
The 1,200—2,400 cycles-per-sec. sound levels did not de- 
crease as much as anticipated for positions near the rear 
of the plane because of inadequate seals around the 
edges of the rear entry way. Similar data would be 
obtained forward of the propeller plane. In general, 
the authors’ experiences indicate that the sound levels 
in Region A decrease at a rate of about 1 db. for each 
4 ft. behind the plane of the propellers and at a slightly 
greater rate ahead. However, standing waves in the 
compartment, air leaks, and thin or untreated fuselage 
surfaces can upset predictions. 


V. PROCEDURE FOR DESIGNING ACOUSTICAL 
INSTALLATION 


Basic Considerations 


In the preceding sections of this paper, the funda- 
mental information to be used in designing an acousti- 
cal installation was presented. In this part a detailed 
procedure will be presented which will indicate the 
proper installation of treatment to use, the portions of 
the fuselage which need greater attention, and the 
minimum weight required to achieve a desired sound 
spectrum. The first step is to attempt to influence the 
basic design of the airplane such that the following 
constructional features are incorporated: 
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(1) The minimum propeller-tip clearance should 
be 20 in. or more. 

(2) Collector ring or other quiet types of exhausts 
should be used. 

(3) Whenever possible, the plane of the propellers 
should be separated from the compartment under 
treatment by means of a bulkhead and should be as far 
from the compartment as possible. 

(4) Windows should have surface densities of not 
less than 1.5 lb. per sq.ft. 

(5) All cracks and openings around windows, doors, 
and hatches should be sealed. 

(6) All accessories should be so mounted that the 
acoustical material can cover as large a portion of the 
fuselage walls as possible. 

(7) All secondary sources of noise, such as genera- 
tors, compressors, etc., should be kept out of the com- 
partment. 


Calculation of the Untreated Noise Spectrum 


To demonstrate the use of the theory, an example for 
a hypothetical transport airplane will be worked out. 
Assume that the flight deck is to be treated and that 
the following information is available: 


Number of blades per propeller a 
Propeller-engine gear ratio 1:2 
Engine r.p.m. at normal cruising 1,800 
Propeller diameter 16 ft. 
Minimum propeller-tip separation 12 in. 


Power per engine at normal cruising 1,000 hp. 
Propeller-tip speed = (1,800 X 16/60)(1/2) = 
754 ft. per sec. 

From Fig. 13, multiplying the abscissa of (a) by */, 
to convert from three to four blades and taking into 
account the tip separation, the level in Region A is 
ound to be 109 db. The octave sound levels in this 
region should lie within +3 db. of this value. The 
datum point just computed for Region A is plotted in 
Fig. 26. 

The construction of the forward part of the airplane 
is shown in Fig. 21. The construction of the bounding 
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Fic. 21. Construction of flight deck of hypothetical transport 
used in calculations. 
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TABLE 2 
1 2 3 + 5 6 7 8 9 10 1l 12 13 14 15 16 17 
S, 
Sq.- 
Item Material Ft In. Sq.Ft. 70.3 ca] 7 Sro.3 Sn Srs 0.3 a as Sars Sa: Sas 
1 Dural 100 0.05 0.7 0.0043 0.00039 3 X 10-%* 0.43 0.039 .0003 0 0 0 0 0 0 
2 Dural 100 0.035 0.5 .0085 .00075 3X 10-5 0.85 0.075 0.0030 0 0 0 0 0 0 
3 Dural 20 0.030 0.43 .02 0.0010 4X 10-5 0.24 0. 0.0008 0 0 0 0 0 0 
4 Dural 100 0.025 0.35 0.018 0.0016 6 X 10-5 1.80 0.160 0.0060 0 0 0 0 0 0 
5 Plexiglas 40 0.25 1.65 .0008 0.00007 2.8xX10-* 0.03 0.003 0.0001 0 0 0 0 0 0 
6 Plate Glass 16 0.25 3.0 0.00023 0.00002 8 x 107 otk ~o 7h 0 0 0 0 0 0 
7 Dural 0.02 0.285 0.026 0.0008 5 X 10-8* 0.52 0.016 .0001 0 0 0 0 0 0 
8 Dural 50 0.02 0.285 0.026 0.0008 5 X 10-46 1.30 0. 0.0002 0 0 0 0 . 0 
T = Total Sr 5.17 0.353 0.0105 
15 20 25 


A = Total Sa 


* See text for explanation. 


surfaces of the flight deck is assumed to be as fol- 
lows: 

Floor: (1) Dural, 0.05%in., braced with '/2-in. wide, 
3-in. deep U braces, spaced 6 in. on centers. Total 
area 109 sq.ft.; unbraced sections of dural, 100 
sq.ft. : 

Plating: (2) Area of 0.035-in. dural below window 
line is 100 sq.ft. (3) Area of 0.030-in. dural above the 
instrument panel is 20 sq.ft. (4) Area of 0.025-in. 
dural above window line is 100 sq.ft. 

Windows: (5) Area of 0.25-in. Plexiglas side and 
overhead windows is 40 sq.ft. (6) Area of 0.25-in 
laminated plate glass windshield is 16 sq.ft. 

Front Bulkhead: (7) Area of 0.02-in. dural between 
flight deck and untreated nose compartment is 20 sq.ft. 
(8) Area of 0.02-in. dural between flight deck and the 
mail compartment is 50 sq.ft. 

Crew Members: (9) Pilot, copilot, navigator, radio 
operator, and flight engineer. 

The first step in the procedure is to determine the 
total transmission (7) and absorption (A) of the 
bounding surfaces. By doing this it will be possible to 
determine the relative amount of treatment which must 
be applied to the various surfaces. The computations 
are made at three frequencies—300, 1,000, and 5,000 
cycles per sec., and the results are shown in Table 2. 
In column one of Table 2 the items are listed as enumer- 
ated above. In columns two and three the basic ma- 
terials and their areas are shown. In columns four and 
five the thicknesses and surface densities are given. 
In columns six to eight the transmission coefficients at 
300, 1,000 and 5,000 cycles per sec. are tabulated.* 
In columns nine to eleven the Sr values are given. In 
columns twelve to 14 the absorption coefficients 
at 300, 1,000, and 5,000 cycles per sec. are tabu- 
lated. Finally, in columns 15 to 17 the Sa products are 
listed. 

The transmission coefficients fer those bounding sur- 
faces that are not outside surfaces of the airplane need 
special handling. At 5,000 cycles per sec. 7; should be 
divided by a factor of 20 if the adjoining compartment 
is well sealed but not acoustically treated. A larger 
factor should be used if the compartment is treated 


(between 20 and 100), and a smaller factor if the com- 
partment has air leaks to the outside. For 7; these 
values can be assumed to be about 10 and 40. Although 
low values of 7; for these surfaces violate the assump- 
tion of uniformity of the boundaries in the compart- 
ment, this difficulty will disappear after treatment is 
added to the other surfaces. 

For Item 1 of this example, a poorly sealed lower 
baggage compartment is assumed. At 300 and 1,000 
cycles per sec. the transmission coefficient is determined 
from Fig. 22 which applies to single impervious septa. 
At 5,000 cycles per sec. the same procedure is used 
except that the transmission coefficient is reduced still 
further by a factor of 5 to take account of the attenu- 
ating effects of the compartment below. For Items 2 
to 6, the transmission coefficients at all three frequencies 
are determined from Fig. 22. For Item 7, the front 
bulkhead, it is assumed that the adjoining compartment 
is untreated but sealed. Hence, 7; for it will be divided 
by a factor of 20, 71 by 4, and 793 will be unchanged. 
For Item 8, the same consideration will hold. If, how- 
ever, Item 8 opened into an acoustically treated com- 
partment, the value of 7; would be divided by 100, 
71 by 20, and 70,3 by 8. 

For Item 9, it is customary practice to assume that 
the equivalent Sa value for a man is between 3 and 5. 
The exact values depend on the type of clothing worn 
and on the frequency, and here the values used will be 
Sao; equal to 3, Sa; equal to 4, and Sas equal to 5 for 
each man in the compartment. 

From the totals shown at the bottom of Table 2, the 
A/T ratios for this example are: 


Ao.3/To.3 = 15/5.17 = 2.9 
= 20/0.353 = 57 
A;/T; = 25/0.0105 = 2,380 


To facilitate computation and to illustrate the toler- 
ances that are applicable to the method, three computa- 
tional charts for determining the spectrum levels rela- 
tive to the level in Region A are given in Figs. 23 to 25. 
The spectrum determined using these charts for this 
airplane with no treatment installed and assuming no 
local sources of noise in the compartment, is shown for 
the “standard position”’ in Fig. 26. 
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Fic. 22. Transmission coefficients, r, for three frequencies vs. ¢, surface density of a single impervious septum (window, sheet of 
or plywood, etc.). 


Planning the Acoustical Lining 


If the design requirements for this airplane required 
ability to converse at 2 ft. at a moderately loud voice 
level, the ‘‘speech interference level’’—i.e., the average 
of the levels in Region S, must be reduced by at least 
15 db. This means a change in the A/T ratio at 1,000 
cycles per sec. from its untreated value of 57 to a value 
of about 2,000. Inspection of Table 2 shows that the 
chief offending surfaces (excluding the floor) are Items 
2, 3,4, and 8. Asa quick check, to see if treating these 
surfaces alone will be adequate, assume that 90 per cent 
of their area can be covered with a material such that 
the resulting structure will have an a of 0.7 and a 7; of 


0.00025. The aS for Items 2, 3, 4, and 8 will become 
170. Hence, the total A; (counting the five men) will 
be 190. Sr, for the untreated 10 per cent of Items 2, 3, 
4, and 8 will be 0.029; for Items 1, 5, 6, and 7, 0.058; 
and for the treated 90 per cent of Items 2, 3, 4, and 8 
it will be 0.061. Hence, the total 7; will be 0.148. 
This gives us an A,/7; = 190/0.0148 = 1,280, which 
falls short of the required 2,000. To accomplish our 
goal, try covering the floor with a carpet, thus reducing 
71 for Item 1 to 0.00025 and increasing a; to 0.2; and 
covering the remaining 10 per cent of Items 2, 3, 4, and 
8. These changes give a new A/T; ratio of 229/0.111 
equal to 2,060, which is adequate. 
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Fic. 24. Computational curves for determining sound levels in 1,200-2,400 cycles per sec. octave band given calculated values of 
A/T at 1,000 cycles per sec. 
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Fic. 25. Computational curves for determining sound levels in 4,800-9,600 cycles per sec. octave band given calculated values of 
; A/T at 5,000 cycles per sec. 


In Fig. 27, which is made up of data in reference 3, it 
is shown that of the four structures, C comes nearest to 
giving the desired value of 7:1. The relation between at- 
tenuation in decibels, as given on the ordinate of Fig. 
27, and 7 at any frequency can be found by comparing 
the right- and left-hand ordinates of Fig. 16. For ex- 
ample, an attenuation of 20 db. corresponds to a 7 of 
0.01. 


Calculation of the ‘‘Treated’’ Spectrum 


Having determined by an estimate that the changes 
in treatment discussed above may be adequate, the 
spectrum for the treated case should be calculated in 
detail. The procedure is the same as before and the 
results are shown in Table 3. The A/T values are: 


. Aos/To.3 = 214/6.87 -= 31.2 
229/0.105 = 2,180 
A;/Ts 318/0.0008 4 X 105 


The new-plotted spectrum is shown in Fig. 28. If 
further reduction in levels in Region § is desired, it is 
seen from Table 3 that 7; for Items 1 and 7 will need to 


-be reduced along with further reduction of 7, for Items 


2,4, and 8. The fact that the Sr; values for the large 
areas are about alike is an indication that efficient use 
has been made of the materials. 
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Fic. 26. Calculated spectrum of noise in untreated hypothetical 
transport. 


VI. MISCELLANEOUS CONSIDERATIONS 


Effect of Acoustically Weak Areas 


A study of Table 3 indicates that a small area of the 
fuselage could be made of lightweight material, say an 
area of 1 sq.ft. of 0.02 dural, without affecting notice- 
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Fic. 27. Typical attenuation curves vs. frequency for four 
varieties of acoustical structures suitable for lining aircraft 
cabins. (See reference 3.) 


ably the fofal transmission T. However, the tempta- 
tion to leave small areas untreated will vanish when it 
is remembered that the procedure given in this paper is 
useful only in estimating the average noise spectrum at 
the ‘‘standard position” in a compartment that has 
fairly uniform acoustical boundaries. Actual measure- 
ments in the vicinity of an acoustically weak area will 
show that sound levels increase markedly as a weakened 
area is approached. This fact can be shown approxi- 
mately by theory. 

As an example, consider the case of a single, thin 
window above the head of the pilot in our hypothetical 
transport. If it is assumed that the cabin is well 
treated so that the average level in Region S at a dis- 
tance from the thin window is 75 db., then by inverse 
square law calculations the contours of sound levels 
in Fig. 29 will be arrived at. 

Near the thin window the sound level will be approxi- 
mately 100 db., one that is almost intolerably high in 
the 1,200-2,400 cycles-per-sec. octave band. At the 
pilot’s head the level will be 85 db., while in the center 
of the compartment the level will decrease to 80 db. 
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Fic. 28. Calculated spectrum of noise in treated hypothetical 
transport. 
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Fic. 29. Calculated noise levels produced by a thin window or 
panel above the pilots’ heads. 


The thickening of this one thin window would reduce 
the level to the 75-db. average in the cabin, which was 
to be expected from the computations. This discussion 
applies equally well to an untreated dural panel or to 
cracks around sliding windows. 


Partial Area Coverage 


It is not always possible to cover 100 per cent of the 
fuselage area with the acoustical material because of 


TABLE 3 
1 2 3 4 5 6 7 9 10 ll 12 13 14 15 16 17 
¢, 
Sq.- ¥. Lbs. per 
Item Material Ft. In. Sq.Ft 70.3 1 Sro.3 Sn Sts a as Sao: Sa Sas 
1 Carpet 100 §=0.08 1.2 0.0018 0.00025 1 X 10-* 0.18 0.025 0.0006 0.1 0.2 0.5 10 20 50 
2 C of Fig. 27 100 3 0.65 0.013 0.00013 4x 1078 1.3 0.01 tan 0.7 0.7 0.9 70 70 90 
3 C of Fig. 27 20 3 0.58 0.017 0.00017 5 x 107% 0.34 0.035 — 0.7 0.7 0.9 14 14 18 
4 C of Fig. 27 100 3 0.51 0.025 0.00025 8 X 107% 2.5 0.025 “oe 0.7 0.7 0.9 70 70 90 
5 Plexiglass 40 0.25 1.65 0.0008 0.00007 2.8x10-§ 0.03 0.003 0.0001 0 0 0 0 0 0 
6 Plate Glass 16 «60.25 3.0 0.00002 8 xX 107 é i 0.00001 0 0 0 0 0 0 
7 Dural 20 0.02 0.285 0.026 0.00080 5 xX 10-6 0.52 0.016 0.0001 0 0 0 0 0 0 
8 ed Fig. 27 50 3 0.44 0.040 0.00040 10 -8* 2.7 0.020 vove 0.7 0.7 0.9 35 20 45 
t = Total Sr 6.87 0.105 0.0008 
A = total Sa 214 229 «6318 


* See text for explanation. 
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TABLE 4 
Total 
Added 
Per Cent Coverage Weight, 
of Items 2, 3, 4, Sr. Sn Sts Saozs Sar Sas Lbs.* 
0 4.19 0.295 0.0100 0 0 0 
25 4.67 0.236 0.0075 47 47 61 7 
50 5.16 0.178 O. 94 122 15 
75 5.65 0.111 0.0025 142 142 182 22 
100 6.14 0.061 0 189 243 


Items 1, 5, 6,7, and 9 0.73 0.044 0.0008 25 40 75 27 


OCTAVE SOUND LEVEL IN DB. 
Flat weighting network 0.0002 dyne/emt ref. level 


* Not including trim cloth and fastenings to cover an area of 290 sq.ft. pa ee ~ caae ieee 
TABLE 5 aa 
Per Cent Coverage of Items 2, 3, 4, 8 0 20 40 60 80 100 
300 128-2 Fic. 30. Variation of sound levels in three frequency bands sta 
6,950 16,400 34,000 78°00 105 asa coverage of cabin side walls with 
apr 
control cables, accessories, etc. To determine how cov- sorption and attenuation play an impor tant part and Fur 
" , careful study should be made of a companion paper con 
ering various percentages of the total area affects the : : 
> - : (reference 3) to learn optimum methods of selecting cori 
levels in Regions B, S, and H, computations were made al 
for the hypothetical transport assuming various per- Te ; oats 
centages of coverage for Items 2, 3, 4, and 8. These Pc eel wh 
items include the surfaces on which Structure C of Fig. ext 
27 was used. A tabulation of the results of the computa- The authors wish to acknowledge the generous help of “ 
# tions is shown in Tables 4 and 5. Plots of the sound personnel of the Army Air Forces, Air Matériel Com- it 
oy levels in the 300-600, 1,200-2,400, and above 4,800 mand, Wright Field; of the Navy Bureau of Aero- § pro 
octaves as a function of percentage coverage are shown nautics; and of the many aircraft and acoustical spe 
in Fig. 30. It is seen that as large a coverage as possible materials companies. In addition, they are indebted wit! 
: is desirable, especially at the higher frequencies, where to various members of the staff of the Electro-Acoustic a 
covering the last 25 per cent is as effective as covering Laboratory at Harvard University for suggestions: and a. 
the preceding 30 to 50 per cent of the area. Both ab- assistance. 2,05 
it u 
400 
REFERENCES inte 
pe 1 Beranek, L. L., Nichols, R. H., Jr., Rudmose, H. W., Sleeper, dra; 
re H. P., Jr., Wallace, R. L., Jr., and Ericson, H. L., Principles of basi 
=. Sound Control in Airplanes, O.S.R.D. Report No. 1543, 388 pp., 
1944 (now out of print). 
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The Jet Airplane Utilizing Boundary 
Layer Air for Propulsion 


A. M. O. SMITH" anp HOWARD E. ROBERTSt 
Douglas Aircraft Company, Inc. 


SUMMARY 


Recent studies in boundary layer theory have shown that the 
stability limit of a laminar boundary layer on a flat plate may be 
expressed as a displacement thickness Reynolds Number. These 
investigations have shown that the natural stability limit value of 
approximately 420! may be increased to 70,000? by the application 
of a weak, continuous removal of the boundary layer by suction. 
Furthermore, the boundary layer displacement thickness with 
continuous suction has been shown never to exceed the value 
corresponding to the upper stability limit and thus is seen to be 
stable to infinite plate chord Reynolds Numbers. Therefore, 
recent theory has indicated the exceedingly important fact that, 
whereas a natural laminar boundary layer is definitely limited in 
extent, this limit can be entirely removed by use of continuous 
boundary layer removal. 

In the light of these results, possibilities of large drag reductions 
are pointed out and a discussion of several of the aerodynamic 
problems of removing the boundary layer is presented. A 
specific study is included comparing a large conventional airplane 
with the boundary layer control jet type, and pointing out the 
salient features of this type. 

The boundary layer control jet airplane considered is found to 
excel the conventional in its pay-load carrying ability up to about 
2,020 miles, whereas a jet airplane with ramming intake can excel 
it up to only 1,310 miles. These ranges are cruised at 430 m.p.h., 
400 m.p.h., and 200 m.p.h., respectively, for the boundary layer 
intake jet, ramming intake jet, and conventional airplane. The 
drag reductions due to the removal of the boundary layer are 
based upon existing test data. 


NOMENCLATURE 


AR = aspect ratio 
Cc = Cov/R,, suction quantity parameter 
Cp = airplane drag coefficient 
Cop = profile drag coefficient; Cpp,, without suction corre- 
sponding to drag measured on balance; Cpp,, 
with suction corresponding to drag determined 
by wake survey; CDp,, with suction correspond- 
ing to drag measured on a balance 
G = skin friction drag coefficient 
CL = airplane lift coefficient; Crimes, maximum lift 
coefficient 
suction quantity coefficient; Q@ = where S 
is characteristic area; Cow, quantity coefficient 
based on wing area; Cg5*, quantity coefficient 
based on displacement thickness of boundary 
layer, = 
= Oswald’s span efficiency factor 
= airplane equivalent parasite area, sq.ft. 
= gravitational constant, ft. per sec.? 
= stream total head, Ibs. per sq.ft.; Hi, total head at 
entrance to slot; H:, total head at exit of slot 


Ce = 


Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
July 18-19, 1946. 

* Design Research Engineer, El Segundo Plant. 

t Aerodynamics Engineer, El Segundo Plant. 
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= airplane wing loading, lbs. per sq.ft. 

power, ft.lbs. per sec. ; 

= pressure, lbs. per sq.ft.; pe, local static pressure at 
edge of boundary layer; p,, free stream static 
pressure; p,, static pressure in duct at station z 

total quantity of flow, cu.ft. per sec.; Q., quantity 
flowing into slot, per running ft. 

velocity head, Ibs. per sq.ft.; g., velocity head in 
slot; go, iocal free stream velocity head; ¢q, air- 
plane velocity head 

wing Reynolds Number; R.z, x Reynolds Number, 
Ux/v; Rs, displacement thickness Reynolds 
Number, U95*/» 

wing area, sq.ft.; Sw, wetted area 

power specific fuel consumption, Ibs. per b.hp.-hr. 


R = 


= 
(S.F.C.)p= 


(S.F.C.)r= thrust specific fuel consumption, lbs. per hr. per Ib. 
of thrust 
7 = thrust, lbs. 
U = velocity, ft. per sec.; U., effective turbojet unit in~ 


let velocity as determined by momentum reac- 
tion; U;, turbojet unit exhaust velocity; U,, 
velocity through slot; U,, duct velocity at station 
z; Up, local velocity at edge of boundary layer; 
Ug, flight velocity 

= local velocity in boundary layer, ft. per sec. 

airplane weight, lbs.; Wo, initial weight; W,, final 

weight 

Wa = engine air consumption, Ibs. per sec. 

x = distance parallel to airflow, ft.; x,, slot width 

y = distance normal to boundary layer, ft. 

2 = distance along duct or slot, ft. 

a = slot loss coefficient (slot loss = aq,) 

B = expansion loss coefficient [expansion loss = (p/2)8 X 
(U, — Us)*) 

= boundary layer displacement thickness, ft. 

n = propeller efficiency 

r = parasite drag reduction factor 

v = kinematic viscosity, sq.ft. per sec. 

p = mass density, slugs per cu.ft. 

o = density ratio 


Subscripts 
i = jet airplane 
? = propeller airplane 


INTRODUCTION 


URBOJET ENGINES may be utilized to propel air- 

craft by at least two different methods. The 
method currently in use treats the engine as an entity 
more or less isolated from the rest of the airplane. This. 
is probably a natural consequence of reciprocating- 
engine installation practices. All installation design. 
effort on current aircraft is directed toward the achieve~ 
ment of maximum thermodynamic efficiency of the 
engine—specifically by striving for the lowest duct loss 
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and highest ram. Great effort is expended in order 
to obtain an installation that produces minimum 
aerodynamic changes in the flow field about the air- 
plane. 

The alternative method capitalizes on the fact that 
the great air consumption of turbojet engines may be 
used to modify and improve the flow over the airplane. 
The ram pressure of this second system, however, may 
be extremely low. The principal difference in the 
objective of these two methods is that whereas in the 
former case a thermodynamically ideal installation is 
obtained, the latter case yields an overall aerodynami- 
cally ideal installation at some expense of thermody- 
namic efficiency. Possibilities of the latter system are 
the subject of this paper. 

The basic thrust equation of a jet engine is the 
familiar relationship 


T = (Wa/g)(Us —*U,) 


It is seen from this equation that reductions in L’, tend 
to increase the thrust by reducing the quantity (W,/g) X 
U,, which is the momentum reaction of the air entering 
the engine. U, will be lower than the flight velocity, 
L’,, if the air is taken from the boundary layers of the 
airplane. At the same time, however, the ram will be 
less, and -this causes the jet velocity and airflow to 
decrease. 

Studies of the use of boundary layer air to improve 
the fuel economy of jet aircraft were begun by the 
authors in 1944. The original analyses were made using 
the simplifying but naive assumptions that boundary 
layer air was a source of low momentum air and that no 
airplane drag changes resulted from its induction. Pre- 
liminary studies showed that in order to satisfy the air 
requirements .of the jet engine, the mean velocities of 
the boundary layer air that is inducted are fairly high: 
on the order of one-half of the free stream velocity. 
However, in spite of these assumptions, the calculations 
showed an improvement of cruising fuel consumption 
per mile of between 5 and 10 per cent over the direct 
tam inlet jet airplane. This gain was not so great as 
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hoped for and hardly enough to justify the engineer- 
ing expense of developing a design along these prin- 
ciples. 

Recently the investigation was expanded to consider 
not only inducting air with low momentum reaction but 
also to consider the possibilities of reduction of skin 
friction drag by proper control of the boundary layer. 
During the interim, results of a theoretical investigation 
in Germany” were made available, which showed that 
large reductions in skin friction drag coefficient were 
possible on a flat plate by using a weak continuous 
boundary layer removal to extend the laminar boundary 
layer stability limit, i.e., delay or prevent transition to 
turbulent boundary layer. By control of the boundary 
layer, together with using air with minimum momen- 
tum, considerable improvements in jet aircraft fuel con- 
sumption or range appear possible. This paper indi- 
cates these improvements, discusses a few of the many 
problems in the design of such an airplane, and presents 
a comparative study of the boundary layer intake jet, 
the ramming intake jet, and the reciprocating-engine 
airplane. The study of the first-mentioned type of air- 
craft is based upon such meager experimental data on 
boundary layer control as are available to date. 


FUNDAMENTALS OF JET AIRPLANE CRUISING 
PERFORMANCE 


A turbojet engine differs from a propeller-driving 
engine in that its fuel consumption per pound of thrust 
is substantially constant. Ram pressure helps its 
specific fuel consumption only by effectively increasing 
its compression ratio. This tends to emphasize the 
fact that turbojet engines designed for use in boundary 
layer aircraft should be designed with considerably 
higher compression ratios than current design practice. 
Since thermal efficiency does not vary rapidly with com- 
pression ratio provided the ratio is reasonably high, 
ram improves fuel consumption only indirectly. Sim- 
ilarly, lack of ram does not greatly increase the specific 
consumption. However, as the basic thrust equation 


TABLE 1 


{1) Instantaneous fuel consumption—miles per Ib. 


{2) Condition for best mileage 


(3) Relation between total drag and parasite drag for maximum 


range—parabolic polar 


(4) Lift coefficient for maximum range—parabolic polar 


(5) Range at constant lift coefficient—miles 


Propeller-Driven 


Airplane Jet-Propelled Airplane 
375n_ Cr (1/W) VG 
SF.C)pWC 
( (S.F.C.)r Co 
Ce _ MAXIMUM VG& waxIMUM 
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shows, U, has a direct effect on thrust, and reductions in 
U, lead to improvement in economy. 


Range and Cruising Considerations 


In order to evaluate the gains made by control of the 
boundary layer, a brief examination of the cruising 
equations for a jet and a conventional airplane will be 
useful. These basic relations are shown in Table 1. 

If the specific fuel consumption is entirely inde- 
pendent of speed or power the condition for best mileage 
is given by line 2, Table 1. These values make no 
assumptions regarding the shape of the airplane polar. 
If the airplane polar can be represented by the common 
approximation, Cp = (f/S) + (C,?/weAR), the lift 
coefficient for best economy is given by line 4. Note 
that it is much lower for a jet airplane than for a pro- 
peller-driven airplane, therefore maximum range occurs 
at much higher speeds. These speeds will become still 
higher as f/Sisreduced. Because the optimum cruising 
C, is lower, the value of the drag coefficient for max- 
imum range becomes (4/3)(f/S) instead of the 2f/S of 
the propeller-driven airplane. Induced drag of the jet 
airplane in the maximum range condition is obviously 
one-third the parasite drag, whereas with a propeller 
airplane it is equal to the parasite drag. This relation 
points out the significant fact that the range of a jet 
airplane is affected much more by its cleanness than is a 
propellered airplane. Similarly induced drag is unim- 
portant to it, so that lower aspect ratios are indicated. 

Integration of the expressions in line 1 between 
initial and final weights produces the range formulas of 
line 5. Provided and 3.F.C.p remain constant, the 
propeller (Bréguet) range formula indicates infinite 
range as the fuel load approaches 100 per cent of the 
gross weight. The jet airplane formula, however, re- 
mains finite. Therefore, before extrapolating the con- 
clusions of range-comparison studies of jet airplanes 
with small percentages of fuel load to airplanes where 
the fuel load is a large fraction of the initial weight, 
further investigation of these two formulas is necessary. 
Fig. 1 presents curves showing the ratio of the range 
with any fraction fuel to the range with 10 per cent 
fuel. All airplane characteristics other than fuel load 
are held constant. The curves show that up to about 
40 per cent fuel there is only a small percentage differ- 
ence in ranges if ranges with small fuel loads are equal. 
As 40 per cent fuel represents a high value, it can be 
concluded that a propeller offers no particular advan- 
tage for flight at the extremely light final weights. 
Moreover, as seen from Fig. 1, if the jet airplane, by 
virtue of the lighter power-plant weights, can carry 40 
per cent fuel as against 37 per cent for the propeller, any 
disadvantage it has in the final phases of the flight is 
nullified. 

The foregoing analysis indicates that if conclusions 
regarding relative cruising economy are reached in 
studying short-range jet airplanes, these same conclu- 
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Fic. 1. Variation of absolute range with fuel load. 


sions will hold for long-range aircraft with high fuel 
loadings. 


Relative Benefits of Drag Reduction to Propeller and Jet 

Airplanes 

The instantaneous miles per pound cruising relation- 
ships given in line 1, Table 1, provide the basic equations 
for use in studying the effects of drag reductions on the 
range of the two types of aircraft. It is assumed that 
the drag reductions would affect only f/S and the aspect 
ratio and other variables would remain constant. De- 
noting f/S = A(f/S)o where (f/S)o represents current 
drag coefficients for each type, the two instantaneous 
cruising relationships may be rewritten as 


mics) af A 
Ibs. jet (S/S) 


miles A, 


Ibs. ) propeller (f/ 


where A; and A; are constants that are characteristic 
of the power plant and all airplane variables except 
parasite drag. 

In current practice the propeller airplane develops 
about twice the maximum miles per pound of fuel as the 


jet airplane. Therefore, 
mileage-propeller A; (f/S) Vids (1) 
mileage-jet Ai 


The ratio of the improvement factors in present practice 
is unity by definition, so 


A2/A1 = = 2 
Substitution of this value into Eq. (1) yields, 


mileage-jet 


If it is assumed that the same percentage improvement 
in drag can be made to both types, Eq. (2) may be 
written, 
mileage-propeller 
mileage-jet 


= (3) 
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Fic. 2. Potential flow into a slot (Cosio¢ = 0.3). 


The above relation shows that, as the drag of both types 
is reduced, the difference between the fuel mileage de- 
creases. If the drag could be reduced to one-sixteenth 
of its present value (an unlikely reduction), the fuel 
mileages would be equal for the two types of aircraft. 

This paper is restricted to a consideration of the use 
of boundary layer control to reduce the drag of jet air- 
craft because the turbojet uses the largest quantities of 
air. Further studies, however, may show that pro- 
peller-driving engines, particularly gas turbines, may 
use sufficient quantities of air to utilize boundary layer 
control efficiently. Nevertheless, reducing the drag of 
the propeller airplane will remain much more difficult 
than for a jet airplane because of the effects of the pro- 
peller on the boundary layers. 


AERODYNAMIC PROBLEMS IN REMOVAL OF THE 
BouNDARY LAYER 


Considerations in the Design of a Boundary Layer Jet 

Airplane for Good Cruising Economy 

It was shown in the previous section that cruising 
gains will be great if considerable drag reduction can be 
made and if care is taken to provide engine air with a 
minimum of momentum. Removing part of the boun- 
dary layer reduces its thickness. Since flow stability is 
primarily a function of the boundary layer displacement 
thickness Reynolds Number, it appears likely, there- 


fore, that boundary layer removal, when properly done, 
can delay transition. Such delay of transition has 
been demonstrated experimentally at R = 25,000,000 in 
flight tests. Furthermore, theoretical studies of con- 
tinuous suction? on a flat plate show much more stable 
profiles. The stability limit, in terms of displacement 
thickness Reynolds Number, is found to vary with the 
suction parameter, C. This limit is shown to vary 
approximately between 420 for the Blasius profile 
(C = 0) and 70,000 for the limiting case of C = ~. 
Moreover, with continuous suction, the boundary layer 
reaches a constant thickness and the flow is stable for all 
choru Reynolds Numbers. 

There is a rather extensive body of experiments on 
boundary layer removal to delay separation but little 
on boundary layer removal to delay transition. The 
only successful experiments in delaying transition have 
been by use of several parallel slots rather than by con- 
tinuous suction through a porous surface. Because of 
this fact and because a practical design using several 
discrete slots is much easier to visualize than a continu- 
ous porous surface, discussion of the intake system will 
be limited to slot systems, with occasional reference to 
the theory of continuous removal. 


Flow Pattern Around a Slot 


Fig. 2 depicts the calculated flow pattern around a 
slot in potential flow. It does not at all represent the 
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flow around such a slot in shear flow but nevertheless 
indicates certain characteristics applicable to the shear 
flow problem. First note that there is a definite closing 
streamline and all fluid within it is taken into the sink. 
Also notice that on the surface there is an accelerated 
flow region favorable to laminar stability on both sides 
of the slot. However, the fact that there is an acceler- 
ated region indicates that the process of removal is quite 
complex and cannot be regarded as a simple superposi- 
tion of a line sink on an undisturbed surface, because it 
immediately changes the surface shear distribution. 

Also shown in this same figure is the pressure distribu- 
tion along the closing streamline. It is seen that it has 
an adverse pressure gradient. In a flow with deficient 
energy as characterized by boundary layer flow, such a 
gradient should be regarded with suspicion since it 
may be able to cause separation, particularly in a 
generally decelerating flow field with the right com- 
bination of boundary layer profile and quantity re- 
moved. This phenomenon appears to have been en- 
countered by the authors in tests of boundary layer 
removal on a wing. It is surprising because it contra- 
dicts the usual belief that boundary layer removal al- 
ways reduces separation. An offset or scoop-type slot 
will reduce any adverse gradient, but not to a large 
extent. 

As will be shown presently, high slot velocities are 
necessary for slot flow stability. The velocity require- 
ment fixes the width of the slot. The height of the 
closing streamline is several times the slot width and, as 
aresult, the closing streamline and its pressures are not 
greatly changed by small changes in slot contours. 

Fig. 2 clearly shows the complexity of the flow field, 
since both large curvature variations and large pressure 
gradients may occur well within the boundary layer it- 
self. 


Stability of Boundary Layer Flow into a Slot 

If the flow into a given slot is increased, the closing 
streamline reaches further into the higher energy layers 
of the local flow field and a possibility of unstable flow 
appears. Consider a simplified system, as indicated in 
Fig. 3, with fluid flowing along a flat plate into an infi- 
nite two-dimensional slot and into a large reservoir, or 
duct. 


LAYER 


Fic. 3. Slot stability diagram. 


The head loss across the slot can be written as aq,. 
Hence, the flow equation across the slot is 


H, = H, — ag, (4) 


where H, and Hz are mean heads defined as stream 
power/Q. If it can be assumed that the kinetic energy 
actually entering the slot is the same as that within the 
closing streamline a slight distance away, where small 
changes in suction rates have negligible effect on this 
reference boundary layer, then the power transported 
into the slot is 


P = So’ */2pu*dy (5) 
The quantity per second is 
Q: = So” udy (6) 
The mean head is P/Q, or 
Hy = So” udy (7) 
If the fluid is incompressible, Eq. (7) may be written 
Hy _ So’ (u/Us)dy 


So” (u/Uo)dy 


In Eq. (4), g; may be written as '/.U,?. With a 
quantity Q, inducted into a slot of width x,, U,; = 
Q,/x, or gs = */2p(Q,/x;)". Q, may be expressed, in 
terms of the boundary layer displacement thickness, 
5*, asQ, = CgU»6*. Then g, becomes 


qa = (9) 


Substitute Eqs. (8) and (9) in Eq. (4) to obtain the 
general flow relation through the slot, 


Hz _ So? 
(=) 


Since fo” udy = Q,, Eq. (10) can be written alternately, 


g So” (u/Ur)dy L Jo 
The flow is stable only if Hz decreases continuously 
as the quantity into the slot increases. Eq. (11) isa 
general statement of the flow relation. The second 
term on the right-hand side is a stabilizing factor be- 
cause it is always negative. The first term varies from 
zero to unity as y varies from zero to infinity. Its 
slope is always positive and thus is destabilizing. 
Therefore, flow from a boundary layer into a slot is in- 
herently unstable, and certain losses are necessary to 
stabilize it. 
A calculation of such flow from a Blasius profile has 
been made assuming 


a/(x,/8*)? = 1 
The results of this calculation are shown in Fig. 4. As 


expected, the head into the slot has an unfavorable 
slope and losses are stabilizing. The plot shown be- 
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Fic. 4. Static stability of flow from a laminar boundary layer 
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comes stable at Cg approximately equal to 0.3 or y/é* = 
1.1. The plot also indicates that stable flow with 
some “‘ram’’ can be expected. On the basis of this 
analysis it appears that slots should be so designed that 
the minimum flow condition is just to the right of the 
neutral stability peak. The type of instability dis- 
cussed in this section has been encountered in high 
aspect-ratio boundary layer bleeds used beneath scoops. 

It is believed that this type of analysis can be used to 
determine the minimum pumping power through a slot. 
The total head of the air leaving the slot appears to 
equal pp plus a fraction of the local g. If fo is con- 
siderably below free stream static, pumping power may 
be high. On the basis of present knowledge it cannot 
be said whether or not the requirement for stability is 
entirely compatible with the type flow best suited to 
delay transition. 


Duct Design and Overall Induction Losses 


Fig. 5 illustrates schematically the overall ducting 
system. The flow problem is unique in that it is the 
problem of flow in a duct with a slot along one side 
through which air is entering. It is assumed that two 
losses are sustained by the air in passing from the 
boundary layer into the duct—one that is proportional 
to slot g, and another that is an expansion loss arising 
when the fluid leaving the slit does not have a velocity 
the same as that in the duct itself. The first loss is the 
stabilizing loss discussed in the previous section. 

Let 


pz = static pressure in duct 
U, = reference slot velocity 


slot loss factor, a = 0 corresponds to no losses 
expansion loss factor, 8 = 1 corresponds to 
sudden expansion 


Then, if the fluid enters the duct essentially parallel 
to its mean stream direction, 


= Hy — (ap/2)U.? — (12a) 


a 
B 


or 


U, = V (2/p)((Hi — + @)] 


The loss in expansion from the slot into the duct is 
(p/2)B(U, — U,)*. The total loss then is 


2 pP [ (2 - 
Ui +58 Vv Uz 

(A, — pz) — ?, 

or, in terms of q,: 
AH (HH, — Pp.) U, — p.) 
Yo Ya(l + a) 


B(U,/U«)? (18) 


U, is readily determined as the integral of the inflow 
from the end of the duct to any point z divided by the 
duct area at z. This equation can be used as the basis 
for determining duct and slot relations. If the duct 
area and removal quantity at any wing station are 
known, a solution for slot width can be made by step- 
by-step integration. 

An arbitrary example was studied for a wing with 2- 
to-1 taper and with a duct whose area was a constant 
percentage of each airfoil section area. It was found 
that the slot should be about twice as wide at the tip as 
at the root because of the static pressure variation along 
the duct. Such an arrangement leads to variable head 
loss into the slot. If the stability criterion is to be 
followed, such variation would be undesirable since a 
constant AH/q,, is desired. 

For a slot at constant chord on a wing of constant 
airfoil section the solution is rather simple. Examine 
Eq. (13); @ and 8 should be constant. If stability is 
to be maintained, H, will also be constant because the 
same fraction of the boundary layer is being removed at 
every point. Then if p, is constant, U,/U. constant 
will product a duct that meets requirements. In such 
a duct p, will vary slightly because of wall friction 
losses along its length. However, this loss is not great. 

In any case, the ideal duct appears to be one whose 
area is proportional to the integral of the flow quantity 
beginning at its extreme tip, and the slot width will then 
be proportional to wing chord. The design criterion 
for ducts with more complicated geometry should al- 


‘ways be to meet the slot stability requirement, which is 


not necessarily constant along the span if the local g or 
boundary layer profile is variable. In brief, the accu- 
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rate design of a slot-duct system is seen to be extremely 
complex. It is likely that an airplane will require 
an extensive tuning up process in a wind tunnel or 
in flight to get the optimum distribution of induced 
air. 


Flying Qualities of Jet Aircraft with Boundary Layer 
Removal 


The anticipated flying qualities of a jet airplane that 
continuously withdraws air from its boundary layers 
must be predicted from scattered and unrelated tests 
of a few aircraft, and principally from component tests. 

The low-speed full-power stability and handling 
qualities may be predicted with more confidence than 
the high-speed characteristics because of extensive flight 
tests conducted on two conventional aircraft by the 
Germans. Although these tests were conducted on 
propeller-driven aircraft in which the blower capacity 
was inadequate, considerable experience and valuable 
data were obtained. 

The behavior of boundary layer controlled aircraft in 
power on stalls leaves much to be desired and points to 
a definite need for an automatic stall-warning device 
coupled with the air-speed indicator and flap position. 
Since strong boundary layer suction leaves little wake, 
no stall warning is given, and in every such airplane a 
violent wing-dropping or undamped rolling oscillation 
about the longitudinal axis has followed the stall. A 
description of a continuously maintained stall empha- 
sizes these adverse characteristics: ‘“‘At slow stalling, 
the plane first tips to the right by about 10 to 15°, 
then to the left with a sudden jump to about 70° lateral 
inclination, then follows a very rapid displacement to 
the right to about 160° lateral inclination. Thus the 
plane was almost on its back, and was again brought 
into its normal position after a half-loop downwards.” 
When not held continuously in the stall, stalling and 
recovery are normal except for a greater loss of alti- 
tude. 


forms and spanwise suction distribution will improve 
the stalling characteristics of suction wings. An inci- 
dental and significant fact in this connection is that 
these adverse stalling characteristics were predicted on 
the basis of wind-tunnel tests. 

Fortunately, no adverse characteristics have been 
observed when the suction fails near the stall. The 
plane merely reduces its angle of attack until flying 
speed is again obtained and no tendency is observed 
either to pitch or roll violently. This is an exceedingly 
important factor in the event of engine failure in take-off 
orlanding. The experience in such failure is likened to 
the experience obtained when ordinary flaps are 
rapidly retracted during flight. 

The take-off and landing characteristics of boundary 
layer control aircraft are expected to be exceptionally 
good. The high maximum lift coefficients that are 
obtained give low take-off speeds and good rates of 
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There is the possibility that improved wing plan 
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climb after take-off. The take-off distances are, 
therefore, considerably reduced. The optimum landing 
flap setting for minimum distance to take-off and climb 
over a 50-ft. obstacle may be considerably greater than 
those normally used for flaps without suction. Section 
profile drags of Cop = 0.005 with take-off suction and 
flaps 60° have been obtained. Minimum landing 
distances may require the use of ‘‘jet thrust destroyer”’ 
since maximum lift coefficients are obtained at maxi- 
mum suction quantities. The ‘thrust destroyer’ is 
also desirable because of the poor acceleration charac- 
teristics of turbojet engines in the event of ‘‘balked”’ 
landings. This implies a requirement for minimum 
time to change the jet to full thrust from a minimum 
thrust condition. 

An interesting and important characteristic of flaps 
with boundary layer control is the fact that a given lift 
coefficient is obtained at a considerably smaller angle of 
attack than that required by normal flaps. This char- 
acteristic is shown in Fig. 6 which shows the longi- 
tudinal angle required to maintain a given air speed, 
with and without suction, for a given flap setting. This 
fact may allow the landing gear length and weight to be 
considerably reduced for boundary layer control air- 
planes. 

At high power and low speed the control character- 
istics of jet aircraft with boundary layer control should 
be exceptional. Section tests of flaps and control sur- 
faces have demonstrated two significant facts: 

(1) In general, the surface is rendered almost twice 
as effective by the removal of the boundary layer at the 
leading edge of the control surface. 

(2) The maximum deflection before loss of effective- 
ness is considerably increased (depending on the quan- 
tity removed and the removal points), upon application 
of boundary layer removal. 

Because the control surface effectiveness of surfaces 
with boundary layer removal approach as a limit the 
theoretical values predicted by airfoil flap theory, it 
can be seen that the per cent chord and required forces 
to produce a given effectiveness may be considerably 
reduced as compared with conventional control sur- 
faces. 
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With the possible exception of stalling character- 
istics, no adverse effects have been experienced with 
boundary layer control on aircraft stability. In the 
cases where adverse stability characteristics have been 
noted, they also exist on the same airplane in the 
absence of suction. 

Through proper use of boundary layer removal, it is 
expected that many stability and control problems 
associated with high flight Mach Numbers may be 
eased. The thinner boundary layers accompanying 
suction may also reduce changes in longitudinal trim 
occurring at flight above the critical Mach Number. 
Nevertheless, it should be realized that at high speeds 
suction may not remove enough boundary layer to 
make appreciable changes in flow patterns. 


Effect of Boundary Layer Removal on Maximum Lift 


The principal emphasis of boundary layer control re- 
search in the past has been on obtaining extremely high 
maximum lift coefficients. Toward this end, many of 
the tests have been successful. However, few tests 
have been made on modern moderately thick airfoils 
with flaps. 

A statistical correlation of Cimaz, data for flapped air- 
foils (6¢ = 45°) of moderate thickness is shown in Fig. 7 
as a function of Cg. The values shown are all obtained 
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Fic. 8. Lift-suction polar. 


with plain flaps with boundary layer removal slots on 
the upper surface near the leading edge of the flap. 
This figure also shows the extremely low profile drag 
coefficients that are obtained at high flap settings with 
suction. 

A typical lift-suction polar is shown in Fig. 8 for an 
NACA 23009 airfoil with a 20 per cent flap.’ These 
data are for an airfoil with an aspect ratio of 5.0. Fig.9 
shows the lift curves from these data corrected to infin- 
ite aspect ratio and indicates the large gains in ACLmez. 
with suction. In this case, at de = 45°, the ACimes. 
is almost doubled, upon the application of suction, com- 
pared to the value without suction. 

A sink superposed on the upper surface of an airfoil 
has two effects: first, a change in the basic potential 
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flow field and, secondly, a change in location of the rear 
stagnation point because of the reduction in boundary 
layer thickness. A sink on the upper surface of an air- 
foil affects the velocity at the trailing edge and, hence, 
the circulation since the Kutta condition must be met. 
Therefore it changes the angle of zero lift. The poten- 
tial flow effects of such a sink are® 


AC, = 2Cg V (1 + 2n)/(1 — 2n) 


where n is the location of the slot aft of the 50 per cent 
chord point in chord fractions. As seen from this ex- 


_ pression, a sink does not affect the slope of the lift 


curve; it only offsets the curve. 

Suction causes the boundary layer near the trailing 
edge to be thinner and to have increased velocities near 
the surface of the airfoil. Then the circulation must be 
increased in order to maintain the Kutta condition at 
the trailing edge. These considerations explain the 
marked changes in trim and lift curve slope shown in 
Figs. 6 and 9. , 


APPLICATION OF BOUNDARY LAYER CONTROL 
To ACTUAL AIRCRAFT 


The previous sections have been devoted to con- 
siderations of problems with individual slots. Princi- 
ples for designing these slots into an actual airplane 
will now be discussed. It must be remembered that air 
is being removed solely to reduce the skin friction or to 
reduce momentum of air entering the engines; hence, 
slot design, size, location, and air distribution always 
should be directed toward obtaining the lowest drag 
and lowest engine air momentum loss. Spacing of 
individual slots cannot be foretold, because of lack of 
data, but probably a considerable number will be re- 
quired to maintain laminar flow. It is likely that gains 
will be made by taking other air from the extremely low 
energy boundary layer profiles on the aft portions. 
Hence, another or several slots may be located in this 
region. Flaps and aileron gaps cause slight separations 
which create sources of dead air. More air should be 
removed where there is a bad local fairing, in order to 
eliminate local separation. 

These principles may be illustrated by a hypothetical 
application to a wing design. Assume that natural 
transition occurs at 20 per cent chord of a 20-ft. chord 
wing at 400 m.p.h. at 30,000 ft. However, the airfoil 
shape is such that pressure gradients are stable to 50 
per cent chord. Then, based upon present knowledge, 
six to twelve slots of 0.02 to 0.05 in. width will be re- 
quired from about 20 per cent to 50 per cent chord to 
maintain laminar flow (see Airplane C—Fig. 13). 
Tunnel tests may show additional gains by locating 
another larger slot at 70 per cent chord and another on 
the flapped surfaces after the joint. Slots will be on 
both upper and lower wing surfaces. 

The fuselage must be treated in the same manner. 
Because of its higher Reynolds Number, it promises to 
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show more gains than wings. Until more is known 
about stability of laminar flow in adverse gradients, care 
should be taken in fuselage design to maintain a stable 
gradient as far aft as possible, say to its midpoint. 
Full peripheral slots appropriately located along the 
fuselage then should maintain laminar flow to transition 
Reynolds Numbers of 100,000,000 or greater. 

Considerable study has been made of the air require- 
ments for stability of boundary layer flow but neither 
theory nor experiment is now adequate to indicate the 
answer. The optimum is probably just enough suction 
to stabilize the boundary layer. More would bring in 
air with greater momentum and greater duct loss, and 
reduce engine efficiency. In case the total engine air 
requirements exceed the requirements for getting min- 
imum drag, it may be desirable to operate some engines 
on direct ramming air. 

The effect of boundary layer control on drag can be 
conveniently expressed in terms of an effective inlet 
velocity ratio. The reaction of a sink in potential flow 
in terms of drag coefficient is ACp = 2Cow. Then if 
Cdp, represents the friction and pressure drag on a 
surface with boundary layer control in operation the 
total drag is: 

Copy = Copa + 


If the drag without suction has a value of Copp, the total 
drag can also be written: 
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Copy = Cop, 2Cow( U./ UVa) 


where U,/U., represents the ratio of effective momen- 
tum of the air entering the sink to its momentum in 
potential flow. Combining the above two expressions 
U,/U., the effectiveness of suction in controlling the 
boundary layer, can be written: 


U./U. = (Copy Cop,)/2Cow 


In this case Cop, and Copp, represent the total drag with 
and without suction such as would be indicated by the 
drag balances in a wind tunnel if the induced air leaves 
the model at right angles to the air stream. U,/U. = 
0 indicates that effectively air has been taken into the 
airplane without momentum reaction. Negative values 
of U./U. indicate that marked reductions in skin fric- 
tion drag have been made. Collected results showing 
the effectiveness of boundary layer control are shown in 
Fig. 10. 

The fact that the effect on U,/U. of removing air 
from the boundary layer is much greater than simply 
inducting air at low velocities can be seen by comparing, 
in Fig. 10, the curve for actual force measurements and 
that calculated assuming that the effective inlet 
velocity is given by the integral 
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Performance Characteristics of a Large Transport Aircraft 
with (A) Reciprocating Engines, (B) Turbojet Engines 
with Direct Ramming Inlets, and (C) Turbojet Engines 
with Boundary Layer Inlets 
Boundary layer air intakes may be applied to jet air- 

craft of any size. The greatest benefits to performance 

however, are believed to be obtained on large aircraft 
with a sufficient number of turbojet engines such that 
efficient cruise operation at full power can be obtained 
by shutting down the required number of units. Large 
aircraft are also well suited because (1) they provide 
more free internal area for ducting, (2) the turbojet 

engines may be submerged in the wing and fuselage, (3) 

increased reliability is obtained in the event of an 

engine failure, and (4) they profit more by delay of 
boundary layer transition than small ones. 

For the purpose of comparison of performance charac- 
teristics, three hypothetical airplanes are assumed and 
shown in Figs. 11, 12, and 13. The dimensional and 
physical data for each of the airplanes accompany the 
figures. It is noted that Airplanes B and C have 
essentially the same characteristics and each has the 
same wing area as Airplane A. The wing aspect ratios, 
however, have been considerably reduced because of the 
fact that these airplanes have much higher cruising 
speeds where lower aspect ratios are permissible. 

The performance comparison is made on the basis of 
all three airplanes having a take-off gross weight of 
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170,000 Ibs. and carrying a cargo of 60,000 Ibs. The 
number of jet engines used for Airplanes B and C was 
selected to give essentially the same static thrust as is 
obtained with the propellers on Airplane A. Reduc- 
tions in wing weight and power-plant installation weight 
were taken into consideration, which resulted in a con- 
siderable reduction of the weight empty of Airplanes B 
and C. This accounts for the fact that a larger amount 
of fuel can be carried in these airplanes than in Air- 


plane A for the same cargo weight and take-off gross - 


weight. 

The parasite drag of Airplanes B and C is also less 
than that for Airplane A because of the removal of the 
engine nacelles, cooling drag, and leakage drag. Air- 
planes A and B each have full-span, double-slotted flaps 
which are unnecessary on Airplane C, which uses only 
simple flaps. The difference in parasite drag for Air- 
planes B and C, therefore, is due to wing flap hinge 
fairings and leakage allowance, which are eliminated in 
Airplane C. Compressibility drag rise was taken into 
account in Airplanes B and C and an increase in the 
critical Mach Number of drag divergence of 0.05 was 
made for the Airplane C as compared with Airplane B. 
This allowance is made on the basis of meager drag data 
at high Mach Numbers.’ 


The take-off calculations were based upon a C, for 
take-off and climb over a 50-ft. obstacle of 90 per cent 
Cimaz. for Airplanes A and B and 85 per cent Cimez. for 
Airplane C. Cimez, for Airplanes A and B was assumed 
to be the same (2.06). The Cimaz. for Airplane C, how- 
ever, varies with air speed and power, i.e., Cg. The 
Cimaz. and Cop variation with Cg shown in Fig. 8 were 
used in the take-off calculations. A graphical solution 
is required in order to determine the C,, for take-off and 
climb. The results of the calculations of take-off dis- 
tance over a 50-ft. obstacle at the design. gross weight of 
170,000 Ibs. are 5,180 ft., 5,980 ft., and 4,070 ft. for 
Airplanes A, B, and C, respectively. The required run- 
way length for Airplane B is therefore 15 per cent 
greater than for Airplane A, whereas it is 21 per cent 
shorter for Airplane C at the design take-off gross 
weight. The superiority of Airplane C over either of 
the other airplanes is largely the result of the increase in 
the take-off and climb C, and an increase in the L/D 
ratio during the climb after take-off. 

The calculations of cruising economy and maximum 
range for each of the three airplanes were based upon an 
average weight equal to the take-off gross weight less 
half of the fuel weight. For the purpose of this analysis 
inlet duct losses for Airplane B were assumed equal to 
25 per cent of the duct inlet velocity head, and for Air- 

plane C an absolute pressure ratio of 0.90 at the entrance 
to the compressor is used. Fig. 10 is used to determine 
the sink reaction on Airplane C. The turbojet engines 
assumed are scaled up from current models and the 
specific fuel consumption figures are representative of 
values currently obtained. 
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Fic. 16. Comparison of maximum range vs. pay load. 


The results of cruising economy calculations and 
miles per gallon vs. true air speed are shown in Fig. 14 
for the three airplanes at several altitudes. Gasoline is 
assumed as the fuel for all engines although Airplanes B 
and C could use kerosene or some other fuel. The 
absolute range corresponding to these economies is 
shown in Fig. 15 for the fuel weight carried under the 
assumed loading condition. It is significant to note 
that whereas the speed for maximum range at 20,000 
ft. for Airplanes B and C is essentially the same, Air- 
plane C has a cruising economy 32 per cent greater than 
that of Airplane B, and compares favorably in range 
with that of Airplane A. 
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Fig. 16 shows the variation of absolute range of the | 


three airplanes with pay load. Airplane C is shown to 
excel A in pay-load carrying ability up to 2,020 miles 
range, whereas Airplane B excels A up to only 1,310 
miles range. It is interesting to note that the inherent 
characteristic of jet aircraft as compared with recipro- 
cating engine aircraft is the extreme flatness of the 
cruising economy vs. speed curves. Consequently, this 
type of aircraft may be operated economically over a 
considerably greater speed range than aircraft with 
reciprocating engines. 

A result of secondary interest is shown in Fig. 14 
which shows the maximum speed of Airplane C to be 
22 m.p.h. greater at 20,000 ft. than Airplane B. Both 
of these aircraft are, of course, considerably faster than 
Airplane A. The cruising economy and range of Air- 
planes B and C may be improved slightly by shutting off 
engines at the lower cruising speeds. The average 
value of Cg for Airplane C in cruise is about 0.002. 
Reference to Fig. 10 shows that the U,/U., correspond- 
ing to this Cg is not so low as could be obtained at lower 
Cg values. 

It is possible to further improve the cruising economy 
by a combination of direct ram and boundary layer 
engines. This is possible because with fewer engines 
operating on the boundary layer air, a lower U,/U.. is 
obtained. The remainder of the thrust required to 
maintain the desired speed would then be made up with 
thrust from direct ram engines. 


CONCLUSIONS 


A visualization:of a jet airplane utilizing boundary 
layer air for propulsion has been presented, brief in- 
vestigations of several of the unique aerodynamic prob- 
lems have been outlined and, finally, an estimation of 
the performance of an airplane designed in accordance 
with the previously discussed requirements is presented. 
This estimation is based strictly on the meager test data 
available and carefully avoids estimation based on 
theory. On the basis of test data alone, however, it has 
been shown that propulsion by use of boundary layer 
air has produced a superior airplane. 

Theory is so incomplete that one is hesitant to use it 
as a basis for prediction. Nevertheless, it is interesting 
to consider whether one may proceed much further in 
drag reduction by use of boundary layer control than 
indicated by available test data. Fig. 17 shows flat 
plate laminar skin friction drag coefficients for various 
degrees of continuous suction as calculated by Schlich- 
ting. The drag coefficient of a good present-day air- 
plane is spotted on the curve. Natural transition is 
indicated by the dotted line. Theoretical stability 
calculations indicate a natural barrier to further laminar 
flow at about this line. Hence, airplanes can only 


progress to it and cannot enter the crosshatched region. 
Buszmann and Miinz*® have calculated the stability of a 
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laminar flow with finite suction at infinite Reynolds 
Number and found it stable at all Reynolds Numbers. 
No other calculations are available although theory is 
now developed. However, it is obvious from examina- 
tion of, boundary layer profiles corresponding to the 
various C values of Fig. 17 that they are much more 
stable than a Blasius profile at the same chord Reynolds 
Number. How far the crosshatched region can be pene- 
trated can only be determined by additional theory and 
test. If it can be penetrated halfway to the Blasius 
line the drag of large airplanes will be reduced to about 
one-fourth of present-day values, and fuel consumption 
per mile will drop equally. It is also possible that tur- 
bine-propeller airplanes can profit similarly if their 
engine air consumption is sufficient for stabilizing the 
boundary layer. Probably a porous surface for con- 
tinuous suction will be found less desirable than a 
series of discrete slots. If so, another unknown is how 
near removal in steps approximates continuous removal. 

Apparently a means for crossing the laminar boundary 
layer stability barrier has been found, and the returns 
from penetrating this hitherto forbidden region are so 
great that extensive research should be initiated to learn 
more about this system, its problems and its possi- 
bilities. 
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Aerodynamic Characteristics of 


Rectangular Wings at Supersonic Speeds’ 


E. ARTHUR BONNEYt 
Applied Physics Laboratory, The Johns Hopkins University 


SUMMARY 


By the combined use of the ‘‘Busemann second order approxi- 
mation” theory for the pressure distribution over a two-dimen- 
sional airfoil with a straight leading edge (normal to the flow) 
and the theory for the loss in lift at the tips of a rectangular air- 
foil in supersonic flow, expressions for lift, drag, and moment coef- 
ficients and center of pressure are derived for the rectangular 
type of wing. Booty 

It is shown that, contrary to subsonic theory, the drag due to 
lift of a rectangular airfoil will decrease with decreasing aspect 
ratio; the lift and moment coefficients ‘will also decrease, and 
the center of pressure will move forward: Increased thickness 
ratio will cause the center of pressure ‘to move forward quite 
markedly for airfoils of finite and infinite aspect ratio. Wave 
form drag will vary as the square of the thickness ratio, and air- 
foils of symmetrical cross section with maximum thickness at the 
mid-chord point will have the least wave form drag for a given 
thickness ratio. 

The analysis is limited to angles of attack below which the 
shock wave detaches from the leading edge of the wing. 


NOMENCLATURE 


Basic Symbols 

= wing semispan 

‘= chord 

= aspect. ratio = span?/area = 2b/c 

= Mach Number 

=VM-1 

= static pressure 

= ratio of tip to root thickness for taper in thickness 

= density 

= ratio of specific heats 

= velocity in ft. per sec. 

=, local speed of sound in ft. per sec. 

= dynamic pressure 

= ‘lift 

= drag 

= wing area = 2bc 

Ci = lift coefficient = L/gS 

Cp = total drag coefficient = D/gS 

Cpr = form drag coefficient or wave drag with no lift 

Cow = wave drag-due-to-lift coefficient . 

CG = skin friction drag coefficient per sq.ft. of wetted area 

Cy, = skin friction wing drag coefficient = 2C, 

Cu = moment coefficient about leading edge = m/cgS 
= center of pressure measured from leading edge 

a = angle of attack 
= semivertex angle of leading edge 
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0 = local angle between any point on the surface of the air- 
foil and the free stream direction 

B = Mach angle = sin~! (1/M) 

K; = coefficient for wing form’ drag 

Ky = coefficient for taper in thickness 

t = wing thickness at root 

ty = wing thickness at tip 

T = thickness ratio = t/c 

Subscripts 

U = upper surface of airfoil ° 

L = lower surface 

F = forward half of double wedge airfoil (except when ap- 
plied to drag coefficient) 

= rear half 

0 = free stream conditions (used only where local and free 
stream conditions are discussed together; otherwise 
free stream conditions have no subscript) 

© = two-dimensional or infinite aspect ratio conditions 


INTRODUCTION 


Two-Dimensional Pressure 


ARIOUS METHODS have been developed for calcu- 

lating and presenting values for the pressure on an 
airfoil in two-dimensional supersonic flow where a 
shock wave has formed and is attached to the leading 
edge. The tabular method of reference 1 and the polar 
method of reference 2 are typical examples. However, 
there is no exact expression for the pressure and these 
methods do not lend themselves to convenient handling 
for purposes of developing expressions for the aerody- 
namic coefficients. However, one method that is ex- 


tremely useful and reasonably accurate is the Busemann 


“second order approximation,”’! which is usually ex- 
pressed as follows: 


= 1 + 2y cosec? + y sec*u cosec? X 
(y + cos? 2)6? (1) 


(@ in rad.). 
Reference to the ‘‘Mach triangle” of Fig. 1 indicates 
the following relations: 
1 1 
sin wo = = VM?—1 °°” 
— 1 


Mo @ 


Also, the expression for dynamic pressure can be 
rewritten in a form more convenient for supersonic 
analysis as follows: 
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M-MACH NUMBER 
Fic. 2. Shock wave detachment angle in two-dimensional flow. 


= po Vo?/ 2 
but Vo = aM and a = V ybo/po; therefore 
go = (y/2)poMo? (3) 


Using these relations, Eq. (1) can be rewritten in a 
form more useful to the aerodynamicist : 


b— po _ Ap_ 2 ) 
do 
Ee + — 2)? 
2(M,? — 1)? 


6? (4) 


Ce (5) 


where C; and C, are functions of the Mach Number 
only. 6 is positive when the direction of the free stream 
is toward the surface and negative when away from the 
surface. 

This method is found to be accurate when compared 
to the exact methods! with the error approaching about 
—2 per cent in lift, drag, and moment at an angle of 
attack of about 60 per cent of the angle at which shock 
wave detachment occurs (called detachment angle). 
Above this, the error increases roughly in a parabolic 
manner until at the detachment angle it may amount 
to —10 to —13 per cent, depending on the geometry of 
the airfoil and the Mach Number. The error in center 


of pressure location (neglecting flow separation) when 
compared to the exact methods is negligible right up to 
the detachment angle. 

The angle of air deflection at which the shock wave 
detaches from the nose of the airfoil and above which 
these theories no longer apply decreases with decreasing 
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Mach Number and increasing leading-edge wedge angle,' 
as shown by Fig. 2. 

This method is limited to wings with straight leading 
edges at right angles to the direction of flow and also to 
Mach Numbers above about 1.3. This theory breaks 
down for Mach Numbers below this value, the exact 
point at which it breaks away being again a function of 
the leading-edge wedge angle. 

From Eq. (4), it is possible to determine the aerody- 
namic characteristics of any cross-sectional’ shape of 
airfoil of infinite aspect ratio. If the loss in pressure 
at the tips for finite aspect ratio wings is also known, 
then the characteristics can be found for rectangular 
airfoils of any aspect ratio. For example, for any wing 
of infinite aspect ratio and symmetrical about its chord 
line, the lift coefficient can be determined simply by 
inspection of Eq. (5) as being 


Cr, = 2Cia = 4a/V M? — 1 (6) 
(@ = +a+Q8, depending on the surface.) 


Tip Loss 

The loss in pressure at the tips due to flow from the 
lower to the upper surface has been determined‘ ® and 
is shown in Fig. 3. This figure represents the pressure 
variation from a value of 0 at the tip to the full two- 
dimensional value at the limit of the Mach angle and 
the curve is described by the equation 


Ap/ Ap, = (2/x) V'tan p’/tan 


This correction was obtained by the use of linear con- 
ical flow theory; however, it is assumed in this 
analysis that pressures in the tip regions are modified 
proportionately to give second order results at the 


c 


T 


4 


TAN ’ 
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Fic. 3. Loss at the tip of a rectangular wing (p = actual 
pressure due to tip loss, p. = corresponding two-dimensional 
pressure). 
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> 


Fic. 4. Lift distribution over a rectangular flat plate in super- 
sonic flow. 
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Fic. 5. Lift distribution over a rectangular flat plate at 
R = 1/V M? — 1. 


Fic. 6. Model of pressure distribution of Fig. 5. 


limit of the Mach cones. For high aspect ratios, 
the lift over a flat plate may be represented diagram- 
matrically as shown in Fig. 4. It can be seen from 
Fig. 3, by comparison with a linear loss in lift from 
the limit of the Mach angle to the tip, that the lift in- 
side the Mach cone will be one-half of the two-dimen- 
sional value for an equal area. 

The effect of interference of the two tip effects with 
each other at low aspect ratios has been determined‘? 
and is shown in Figs. 5 and 6. 

It is important to note that the pressure at any point 
in the region of combined tip effects is equal to the sum 
of the pressures as affected by each tip minus the two- 
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dimensional value. Because of this effect, the expres- 
sions developed on the basis of lift distribution at high 
aspect ratios where the tip effects do not cross on the 
wing will also be correct for lower aspect ratios where 
they do cross. This will be true until the aspect ratio 
becomes so low that the limit of the tip loss from one 
tip just reaches the intersection of the trailing edge and 
the opposite tip as represented by the last two figures 
above. This limiting aspect ratio will be 


R =1/VM?—1 


For extremely low aspect ratios, the lift force on the 
“tail” of the wing (aft of c = 2b VM? — 1) will appar- 
ently alternate between relatively small negative and 
positive values; however, because of lack of exact 
quantitative information, the force in this region has 
been assumed to be negligible. ; 

With conditions known at every point on the airfoil, 
it is now possible to develop expressions for lift, drag, 
moment, and center of pressure for any type of airfoil. 


ASSUMPTIONS 


It was necessary to make certain simplifying assump- 
tions to keep the various expressions reasonably simple 
and yet accurate. They are as follows: 

(1) The possibility of secondary tip effects originat- 
ing at the point of maximum thickness of a double- 
wedge airfoil, for example, was not considered. 

(2) The loss in pressure at the tips actually extends 
to the limit of the /ocal Mach angle, which changes on 
each surface with angle of attack. However, the Mach 
angles will be higher on the lower surface and lower on 
the upper surface, thereby offsetting each other to a 
great extent. Therefore, the free stream Mach angle 
is used in the derivations. 

(3) Consideration of the phenomena of separation 
was of necessity omitted because of the lack of informa- 
tion on this effect in supersonic flow. This factor can 
cause the center of pressure expressions obtained herein 


to be somewhat in error, particularly at high angles of |. 


attack. (The reader is referred to reference 3 for a 
study of the effects of separation.) 

(4) Throughout the analysis the assumption is 
made that sin a = a and cos a = 1. 


DERIVATION OF AERODYNAMIC COEFFICIENTS—REC- 
TANGULAR WINGS 


The symmetrical double-wedge section will be used 
in the following illustrative derivations and the final ex- 
pressions generalized for all symmetrical airfoil shapes 
by a method also suggested by Busemann. Only the 
lift and drag coefficients and center of pressure of an 
airfoil of finite aspect ratio will be derived as examples 
of the method involved. 

Fig. 7 represents the cross section and planform of a 
symmetrical (about chord line and midpoint) double- 
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RECTANGULAR WINGS AT SUPERSONIC SPEEDS 


In the region unaffected by tip loss the flow is two- 


ves dimensional, having the value 
Ap/q = Cio + 
® rO\r ~ /© 
—@ 
©) A=-atB 
tas 7 s= atB e-6 


Therefore the two-dimensional pressure difference over 
the front and rear halves of the airfoil will be 


Fic. 7. Cross section and planform of rectangular double-wedge 
airfoil showing areas of tip loss. 


wedge airfoil and shows the limits of the regions of tip 
loss, which are defined by the Mach angle measured 
from the tip. The surfaces and areas are numbered and 


(8) 
lettered for clarity in understanding the derivation to Ap Ap Ap 
follow. —) =({—] = 2Ga — 4C.a8 
q 
Lift Coefficient 


Remembering that the average lift in region C is equal to half of that in region A (for a like area) and likewise 


the lift in region D is half that of B, then 
ifs Cc ac\e 1 /3c 


(2Cia + 4Cra8) | (2 - 


This expression applies to a double-wedge airfoil only. 
However, it has been pointed out’ that the expression 
can be made general for all symmetrical airfoils by sub- 
stituting the parameter of cross-sectional area divided 
by the chord squared in place of the half-wedge angle. 
Therefore, for a double-wedge airfoil, 


Acs/c? = tc/2c? = t/2c = B/2 
therefore, 
B = 2Ags/c? = 2A’ (11) 
Values of A’ for various types of airfoils are noted in 
Table 1. Therefore, if C; = 2/B and C; = 2C./G, 
C, = (4a/B)[1 — (1/2RB)(1 — GA’)] (12) 


Inasmuch as the lift coefficient for infinite aspect 
tatio is C, = 4a/B, then the aspect-thickness ratio cor- 
tection to two-dimensional lift (and drag due to lift) is 


Ci/Cr_ = 1 — (1/2RB)(1— GA’) (13) 


It will be noted that this derivation is for high aspect 
tatios where the tip effects from opposite tips do not 
cross on the wing. As has been mentioned, however, 
the expression will also hold down to a value of aspect 


L/ 
2bc 
(9) 
(where B = V M? — 1). 
Cr = 2Cia E (10) 


ratio = 1/V M? — 1. Below this value the correction 
for aspect ratio will be (assuming lift to be negligible 
aft of the chord point c = 26V M? — 1) 


Ci/Cr. = (RB/2) (1 + GA’) (14) 


It should also be noted that the factor inside the second 
parentheses in these expressions represents a correction 
for thickness ratio, while the main body of the expression 
is the aspect ratio correction for an infinitely thin flat 


plate. The flat plate correction is plotted in Fig. 8, 
R:2.0, 
NESS FLAT PLATE 

Cy 

het | 
4 
2 

RB= 


Fic. 8. Aspect ratio correction to lift and drag due to lift. 
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10 LS 20 25 3.0 35 


Fic.9. Lift curve slope for various values of aspect ratio of a flat 
plate. 


with a 10 per cent double-wedge airfoil of aspect 
ratio = 2 added to show the effect of thickness. 

This curve is replotted in Fig. 9 to show lift curve 
slope versus Mach Number for various aspect ratios of a 
flat plate. The error incurred by using this curve for 
airfoils of finite thickness will be negligible for practical 
values of thickness ratio. 

The factor 


_ + (M? — 2)? 
CG 2(M? 1)” 
is plotted in Fig. 10. 


It should also be noted that the expressions for lift 
coefficient derived above are actually normal force coef- 


C; = 


ficients, inasmuch as the average pressure acts normal © 


to the chord line of the airfoil. However, inasmuch as 
C, = Cy cos a, this effect can be neglected until the 
angle becomes relatively large. 


Drag Coefficient 


The total drag coefficient is equal to the sum of the 
coefficients of wave drag at zero lift (hereinafter re- 
ferred to as form drag) plus the wave drag due to lift 
(wave drag) plus the skin friction. 

By analogy with the derivation for Eq. (6), it is ob- 
vious that the form drag coefficient will be 


Co, = (48/B) sin B = 462/B (15) 


(It is to be noted that this type of drag is unique to 
supersonic flow only, there being no corresponding term 
in subsonic flow—i.e., there is no form drag in a perfect 
fluid subsonically.) 

The drag due to lift will be the dragwise component 
of the normal force coefficient or 


C. 
Coy = sin a = (16) 
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Fic. 10. Mach Number correction factor C;. 


TABLE 1 
Summary of Airfoil Characteristics 
[ type | FLAT | FLAT PLATE | AIRFOL | AIRFOIL 
| R | FINITE INF INIPE | 
4a 4 | 
or A, OL IN RADIANS 
TYPE OF AIRFOIL A 
DOvBLE 4 
533 
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Adding the effect of skin friction, the total drag coef- 


ficient becomes 


48? Cr 
17 
Co + Cp; (17) 
where 8 in the general equation is the effective semi- 
wedge angle of the leading edge. The expression may 
also be written generally in a more usable form as fol- 
lows: 


4a? CL 


Bo B 
Values of K, for various types of airfoils will be found in 
Table 1. 
For the case of a wing tapered in thickness, a correc- 
tion factor is applied to the form drag coefficient: 


K, 


where 7 is the thickness ratio at the root of the wing. 


Co; (18) 
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Fic. 11. Drag correction for taper in thickness. 


The factor K, for wings with single taper is given in 
Fig. 11. 

The equation for the upper curve is K, = 1/2 (1 + A) 
and for the lower curve is K, = 1/3(1 +A + A?2). 

It is obvious from Eq. (17) that an airfoil section that 
is symmetrical about its chord line and midpoint is the 
optimum design, because, for a given thickness ratio 
orstrength, the effective value of 8 will be a minimum for 
this shape. (Although moving the point of maximum 
thickness aft of the midpoint decreases the actual lead- 
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ing edge wedge angle, the effect of the increased trailing 


edge angle is to increase the effective 8.) 


Center of Pressure 


The center of pressure is obtained by summating the 
products of the distances to the center of pressure of 
each region as shown in Fig. 7, by the total pressure 
in that region for the entire airfoil, and dividing by the 


total pressure or lift force. 


The distance to the center of pressure of each area is 


c 6L4b—(c/B) | 12RB—6 
cp.\ 1,1 6b — (5c/B) mM 9RB — 14 
(c.p./c)e = (2/3) (1/2) = 1/3 

(c.p./c)n = 1/2 + (5/9)(1/2) = 7/9 


The area of each region is 


= (2RB — 1) (ses) 


Therefore the center of pressure will be (omitting factor bc/2RB from top and bottom): 


12RB—6 


(20) 


(21) 


Cp. 


(2Cia + 4Cza8)[2RB — 1 + (1/2) (1)] + Gia — 4Cra8)[2RB — 3 + (1/2) (8)] 


c.p./c = [RB — 2/3 — C,A’(RB — 1)]/(2RB — 1 + GA’) 


This expression is likewise general for all symmetrical airfoil sections when expressed as a function of A’ rather 


than 


The variation in center of pressure with the param- 
eter RB is shown in Fig. 12 for a flat plate and for 10 
per cent double-wedge and biconvex airfoils having an 
aspect ratio of 2 as practical examples of the trends to 
be expected. 


It should be noted that this expression makes no al.-. 


lowance for separation effects, which will cause the 
center of pressure to move forward with increased angle 
of attack. Because of the lack of information on the 
magnitude of this effect in supersonic flow, it was not 
possible to incorporate the effect into the expression. 
This theory should, however, aid in determining the 
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Fic. 12. Center of pressure location. 
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magnitude of separation effects from experimental 
data. 

It is also possible that in the actual flow over a 
double-wedge airfoil the boundary layer profile may 
have a curved shape immediately aft of the point of 
maximum thickness. The effective airfoil shape may 
therefore be a continuous curve to the new angle rather 
than an abrupt change. This would decrease the rate 
of expansion on the rear half of the airfoil and cause the 
pressure to be more positive in the region just aft of the 
point of maximum thickness. Because of symmetry, 
the lift force would probably not be affected too 
greatly; the effect on drag would also be small; how- 
ever, the effect on center of pressure could be quite 
marked. 

The aerodynamic expressions for flat plates and air- 
foils both of infinite and finite aspect ratio are listed in 
Table 1, together with values of the various constants 
used throughout this analysis. . 


CONDITIONS 


The previously derived expressions have been used to 
determine optimum conditions of aspect ratio, angle of 
attack, lift-drag ratio, etc., for any given conditions of 
allowable stress, airfoil cross-sectional shape, and Mach 
Number, as well as method of supporting the wing— 
i.e., by the entire base or on a shaft through the body. 
However, this work is omitted here in the interests of 
brevity. It develops that the maximum lift-drag ratio 
that can be expected from supersonic wings of these 
types decreases with increasing Mach Number and is 
generally between one-half and one-quarter of the values 
obtained for practical subsonic wings. 

Because of the excessive weight of solid wings, it is 
probable that some type of monocoque construction 
will probably be the ultimate in design. This would 
mean an increase in the thickness ratio and form drag 
coefficient; however, the decreased wing area required 
as a result of the decreased weight would undoubtedly 
offset this effect and cause the overall drag to be less. 


CONCLUSIONS 


Detailed conclusions to be drawn from this analysis 
are as follows: 

(1) Contrary to subsonic theory, the drag due to lift 
of a rectangular airfoil in supersonic flow will decrease 
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with decreased aspect ratio. The lift and moment coef. 
ficients will also decrease, and the center of pressure 
will move forward with decrease in aspect ratio. 
(2) The center of pressure of a flat plate of infinite 
aspect ratio is at the mid-chord point. 
(3) Increasing the thickness ratio will increase the 
lift and wave drag (due to lift) coefficients of airfoils of 
finite span slightly but will decrease the moment coef- 
ficient (about the leading edge) and cause the center of 
pressure to move forward quite markedly for airfoils 
of finite and infinite aspect ratio. 
(4) Airfoils having the same cross-sectional area will 
have the same center of pressure location. 
(5) Increasing the Mach Number will either increase 
or decrease the coefficients depending on the magnitude 
of the Mach Number being considered and the aspect 
ratio. 
(6) Form drag will vary as the square of the thick. 
ness ratio for a given cross-sectional shape. 
(7) Airfoils of symmetrical cross section with maxi- 
mum thickness at the mid-chord point will have the 
least drag for a given thickness ratio or allowable stress. 
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The Limiting Effect of Centripetal 
Acceleration on Man’s Ability to Move’ 


CHARLES F. CODE,+ EARL H. WOOD,{ ann EDWARD H. LAMBERT, 
Mayo Clinic and Mayo Foundation 


ABSTRACT 


Instances have been recorded in which fliers were unable to bail 
out of spinning aircraft because of the centrifugal force gener- 
ated by the spin. A study designed to give an estimate of the 
restrictions placed on man’s locomotive ability by exposure to 
radial g (centrifugal force) was made on the human centrifuge. 
Five subjects were studied. The average time required to don a 
parachute was increased from 17 sec. at 1 g to 1 min. and 15 sec. 
at 3 g, it being uniformly agreed by the subjects tested that they 
could not don the parachutes at accelerations slightly above 3 g. 
The ability of the subjects to move against the force or at right 
angles to the force was seriously restricted when the magnitude 
of the force was 2 to 3 g. Progress against the force became 
impossible at 3 gz. Moving the body at right angles to the force 
became impossible in the neighborhood of 4 g. The study demon- 
strated that under conditions at which accelerations of 3 g or more 
develop, fliers will need help if they are to escape from their air- 
craft. Emphasis is therefore given to the need, recognized for 
some time by others, of further consideration and development of 
devices designed to assist the flier when escaping from his aircraft. 


INTRODUCTION 


r HAS RECENTLY BEEN EMPHASIZED that escape from 
spinning aircraft may be seriously hampered by the 
forces resulting from the spinning motion of the aircraft.’ 
Exact data on the amount of force arising in different 
types of aircraft during various-spins are not available, 
but it has been pointed out that one of the sensations 
experienced by aviators under such circumstances is 
that an “‘invisible force’’ of great magnitude so presses 
down on them that, when they attempt to escape, their 
bodies are practically immobilized and only their arms 
and hands may be freely moved. 

A spinning aircraft is not unlike a spinning human 
centrifuge. The forces that develop in the two instances 
are similar. The major component of these forces re- 
sults from the centripetal acceleration generated by the 
spinning motion of the plane or the centrifuge. The 
reactive force to this acceleration is the ‘‘invisible force” 
to which fliers have ascribed their difficulties. 

On the human centrifuge the plane of rotation and 
the orientation of the fixed parts of the centrifuge in 
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relation to the force remain constant. In spinning air- 
craft, particularly during violent spins arising from 
structural damage, the orientation of the aircraft and 
the plane of rotation may change rapidly, resulting in 
sudden alterations in the direction of the force with the 
result that the contents of the plane, the fliers included, 
may be thrown about violently within the ship. These 
violent changes in the direction of the force cannot be 
conveniently duplicated on the human centrifuge. 
Exact studies under accurately known and controlled 
conditions, free from the hazards of flight during spins 
can, however, be made on the centrifuge. 

While the exact accelerations generated in spinning 
aircraft are not known, it is known that the flier may be 
pinned down by the force so that he is unable to move 
his body. This study was undertaken with the aim of 
obtaining some measure of the restrictions placed on 
man’s ability to move and carry out certain maneuvers 
by centripetal accelerations of various magnitudes and 
also to determine the acceleration at which effective 
movement of the body is no longer possible. 


EXPERIMENTAL PROCEDURE 


The experiments were carried out on the platform or 
open end of the human centrifuge in the Mayo Aero- 
Medical Unit.? (Fig. 1.) At this end of the centrifuge 
there is an open space approximately 7 by 7 by 7 ft. 
(2.1 by 2.1 by 2.1 m.) which affords ample room for 
tests that involve walking, crawling, etc. The sides and 


Fic. 1. The open or platform end of the human centrifuge 
caged in, reinforced and covered with gymnasium mats so that 
subjects could not walk, crawl, or fall off the centrifuge or injure 
themselves if they fell during tests. 
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CONTROL 


START MIDWAY 


DURING CENTRIPETAL ACCELERATION 
START FINIS AND EXHAUSTION 


3.57 


Fic. 2. Crawling across the end of the centrifuge during centripetal acceleration—moving at right angles to the direction of force. 
In this and subsequent illustrations all pictures are single frames taken from motion pictures. Control observations were made at 1 g 
gravity. Notice clock which gives time in seconds, 
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top of this end of the centrifuge were covered with flat 
expanded metal mesh sheeting so that it was impossible 
for a man to walk, crawl, or fall off the centrifuge during 
tests. The end of the centrifuge was reinforced so that a 
falling man could not break through. 

Tests of the strength of the caging and reinforcement 
were first made using a dummy weighing 200 lbs. (90.7 
kg.). The dummy, supported from the center, was cut 
loose 2 to 4 ft. (61 to 122 cm.) from the end of the centri- 
fuge during accelerations of 2 to 3 g. The procedure 
simulated a man slipping and falling while trying to 
make headway against radial g. In initial tests the 
dummy tore the caging loose and it was necessary to 
reinforce its attachments before proceeding. The 
experience served as a warning of the danger to which 
subjects may be exposed in such tests if the magnitude 
and the effects of these forces are underestimated. As 
some protection to the subjects, the end, and at times 
the sides, of the centrifuge were covered by standard 
gymnasium mats 2 in. (5 cm.) thick (Fig. 1). 

To aid in standardization and simplification of ex- 
pression the only component of acceleration recorded 
in these experiments was the centripetal acceleration 
acting in a direction perpendicular to the end of the 
centrifuge. This acceleration is referred to in this 
report as radial g. It should be pointed out that while 
radial g is the major component of the forces to which 
the subject is exposed on the rotating centrifuge, the 
pull of the earth’s gravitational field (1 g, gravity) is 
also acting. At centripetal accelerations of 2 g and more 
the contribution made by the earth’s field of gravity to 
the resultant is small and for this reason has been 
neglected. Ihe maximal radial g of each run was deter- 
mined by calculation from the r.p.m. of the centrifuge 
and was recorded on a Kollsman accelerometer at the 
end of the centrifuge—a distance of 17*/; ft. (5.4 m.) 
from the center of rotation. 

In order to escape from a spinning aircraft the flier 
may be required to move in any direction with respect 
to the direction of the force. For simplification on the 
centrifuge four types of movement with respect to the 
force were tested: 

(I) Movement in the same direction as the radial 
force (moving with the radial g). No special tests were 
designated for movement in this direction. Some ob- 
servations were made, incidentally, while carrying out 
the other tests. 

(II) Movement at right angles to the direction of 
the force (moving across the radial g). Two tests of 
movement in this direction were made: 

(a) In the first test the subjects were required to 
crawl from a sitting position at one side of the end of 
the centrifuge to a sitting position at the other side 
of the end of the centrifuge, covering a distance of 
7'/. ft. (2.3 m.) (Fig. 2). The time required to accom- 
plish this test at 1 g gravity (three or more trials) 
was compared with the time required at different 
levels of radial g. The results were expressed as a 


percentage increase over the time required at 1 g 

gravity. 

(b) The second test was similar except that the 
subject had to round a barrier placed midway be- 
tween the two sides of the end of the centrifuge (Fig. 
3). The barrier projected 22 in. (56 cm.) from the 
back of the centrifuge. 

(III) Movements against radial g: The direction 
of movement under these circumstances was toward the 
center of the centrifuge. The effect of radial g on the 
subject’s ability to walk, crawl, climb a rope or a ladder, 
or to rise from a bomber seat toward the center of the 
centrifuge was tested. 

(IV) In the fourth test the subject was required to 
carry out the comparatively complicated maneuver of 
donning a standard back parachute. The effect of 
radial g on the time required to complete this test was 
determined. 

Five subjects, all normal men, were studied. Not all 
of the subjects carried out all of the tests. During all 
tests the subjects wore the standard Army flying cover- 
all; in some they also wore electrically heated suits. 

Motion pictures were made during all tests. These 
were subsequently studied, and the results noted at the 
time of the test were checked and rechecked in the mo- 
tion pictures. The film illustrates the difficulties and 
limitations of movement under radial g. During all 
tests the camera was mounted at the center of the 
centrifuge. 


RESULTS 


(1) Movement in the Same Direction as Radial g.— 
Movement in this direction is easy but it may be hazard- 
ous. The first example of its danger was seen when the 
200-Ib. (90.7-kg.) dummy, cut loose 3 ft. (91 cm.) from 
the end of the centrifuge at 2 g, sprung some of the 
caging. When the subjects slipped, stumbled, or fell at 
1 radial g, they hit the end of the centrifuge with a force 
approximately equal to that pertaining had they fallen 
the same distance to the earth. At 2 g the hazards were 
increased, the body weight of the subjects being doubled 
[a 200-Ib. (90.7-kg.) subject now weighing 400 Ibs. 
(181.4 kg.)]; at 4 g, movement in this direction was 
made additionally dangerotis by the fourfold increase in 
body weight. Falling just a few feet under these 
conditions resulted in exceedingly high impact forces 
when the subjects hit the side or end of the centrifuge. 
The same circumstances can and must at times pertain 
in spinning aircraft. These impact forces were not 
measured in this study and movements in this direction 
were not observed at accelerations greater than 2 to 3 
g. Falling or dropping a short distance [1 ft. (30 cm.)] 
at 2 g was sufficient to knock the wind out of subjects. 
It was evident that serious damage could easily be done 
by short falls at accelerations above these values. 

(2) Movement at Right Angles to Radial g.—(a) 
Crawling across the end of the centrifuge, a distance of 
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Fic. 2. Crawling across the end of the centrifuge during centripetal acceleration—moving at right angles to the direction of force. 


In this and subsequent illustrations all pictures are sin 
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top of this end of the centrifuge were covered with flat 
expanded metal mesh sheeting so that it was impossible 
for a man to walk, crawl, or fall off the centrifuge during 
tests. The end of the centrifuge was reinforced so that a 
falling man could not break through. 

Tests of the strength of the caging and reinforcement 
were first made using a dummy weighing 200 lbs. (90.7 
kg.). The dummy, supported from the center, was cut 
loose 2 to 4 ft. (61 to 122 cm.) from the end of the centri- 
fuge during accelerations of 2 to 3 g. The procedure 
simulated a man slipping and falling while trying to 
make headway against radial g. In initial tests the 
dummy tore the caging loose and it was necessary to 
reinforce its attachments before proceeding. The 
experience served as a warning of the danger to which 
subjects may be exposed in such tests if the magnitude 
and the effects of these forces are underestimated. As 
some protection to the subjects, the end, and at times 
the sides, of the centrifuge were covered by standard 
gymnasium mats 2 in. (5 cm.) thick (Fig. 1). 

To aid in standardization and simplification of ex- 
pression the only component of acceleration recorded 
in these experiments was the centripetal acceleration 
acting in a direction perpendicular to the end of the 
centrifuge. This acceleration is referred to in this 
report as radial g. It should be pointed out that while 
radial g is the major component of the forces to which 
the subject is exposed on the rotating centrifuge, the 
pull of the earth’s gravitational field (1 g, gravity) is 
also acting. At centripetal accelerations of 2 g and more 
the contribution made by the earth’s field of gravity to 
the resultant is small and for this reason has been 
neglected. The maximal radial g of each run was deter- 
mined by calculation from the r.p.m. of the centrifuge 
and was recorded on a Kollsman accelerometer at the 
end of the centrifuge—a distance of 174/s ft. (5.4 m.) 
from the center of rotation. 

In order to escape from a spinning aircraft the flier 
may be required to move in any direction with respect 
to the direction of the force. For simplification on the 
centrifuge four types of movement with respect to the 
force were tested: 

(I) Movement in the same direction as the radial 
force (moving with the radial g). No special tests were 
designated for movement in this direction. Some ob- 
servations were made, incidentally, while carrying out 
the other tests. 

(II) Movement at right angles to the direction of 
the force (moving across the radial g). Two tests of 
movement in this direction were made: 

(a) In the first test the subjects were required to 
crawl from a sitting position at one side of the end of 
the centrifuge to a sitting position at the other side 
of the end of the centrifuge, covering a distance of 
7'/. ft. (2.3 m.) (Fig. 2). The time required to accom- 
plish this test at 1 g gravity (three or more trials) 
was compared with the time required at different 
levels of radial gz. The results were expressed as a 
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percentage increase over the time required at 1 g 

gravity. 

(b) The second test was similar except that the 
subject had to round a barrier placed midway be- 
tween the two sides of the end of the centrifuge (Fig. 
3). The barrier projected 22 in. (56 cm.) from the 
back of the centrifuge. 

(III) Movements against radial g: The direction 
of movement under these circumstances was toward the 
center of the centrifuge. The effect of radial g on the 
subject’s ability to walk, crawl, climb a rope or a ladder, 
or to rise from a bomber seat toward the center of the 
centrifuge was tested. ‘ 

(IV) In the fourth test the subject was required to 
carry out the comparatively complicated maneuver of 
donning a standard back parachute. The effect of 
radial g on the time required to complete this test was 
determined. 

Five subjects, all normal men, were studied. Not all 
of the subjects carried out all of the tests. During all 
tests the subjects wore the standard Army flying cover- 
all; in some they also wore electrically heated suits. 

Motion pictures were made during all tests. These 
were subsequently studied, and the results noted at the 
time of the test were checked and rechecked in the mo- 
tion pictures. The film illustrates the difficulties and 
limitations of movement under radial g. During all 
tests the camera was mounted at the center of the 
centrifuge. 


RESULTS 


(1) Movement in the Same Direction as Radial g.— 
Movement in this direction is easy but it may be hazard- 
ous. The first example of its danger was seen when the 
200-lb. (90.7-kg.) dummy, cut loose 3 ft. (91 cm.) from 
the end of the centrifuge at 2 g, sprung some of the 
caging. When the subjects slipped, stumbled, or fell at 
1 radial g, they hit the end of the centrifuge with a force 
approximately equal to that pertaining had they fallen 
the same distance to the earth. At 2 g the hazards were 
increased, the body weight of the subjects being doubled 
[a 200-lb. (90.7-kg.) subject now weighing 400 Ibs. 
(181.4 kg.)]; at 4 g, movement in this direction was 
made additionally dangerotis by the fourfold increase in 
body weight. Falling just a few feet under these 
conditions resulted in exceedingly high impact forces 
when the subjects hit the side or end of the centrifuge. 
The same circumstances can and must at times pertain 
in spinning aircraft. These impact forces were not 
measured in this study and movements in this direction 
were not observed at accelerations greater than 2 to 3 
g. Falling or dropping a short distance [1 ft. (30 cm.)] 
at 2 g was sufficient to knock the wind out of subjects. 
It was evident that serious damage could easily be done 
by short falls at accelerations above these values. 

(2) Movement at Right Angles to Radial g.—(a) 
Crawling across the end of the centrifuge, a distance of 
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71/2 ft. (2.3 m.). Five subjects were tested. On the 
average the time required for this maneuver was in- 
creased two and one-fourth times by 1 g radial, five times 
by 2 g radial, and nine times by 3 g radial (Table 1). 
Three of the five subjects tested could not accomplish 
the task at 4 g (Fig. 2). At or slightly above 4 g all 
of the subjects were immobilized. They could move 
their arms and legs, and these they pushed back and 
forth against the end of the centrifuge in an apparently 
aimless fashion. Their bodies were pinned to the end 
of the centrifuge with such force that they simply could 
not make any progress. While the gymnasium mat at 
the end of the centrifuge provided good “footing,” it 
also offered considerable friction to body movement. 
With suitable hand and foot grips some progress at 4 g 
might have been possible had a smoother, harder sur- 
face been provided. At 4 g, while forward or backward 
crawling was impossible, some of the subjects were 
able to roll themselves up the back of the centrifuge 
(Fig. 4). This motion might be of value in reaching an 
escape hatch. 

(b) Rounding a barrier. The motion required to get 
past the barrier was not unlike that which would be 
needed to get by a bulkhead in an aircraft (Fig. 3). 
The time required to round the barrier was increased 
two and one-fourth times by 1 g radial, six times by 2 g 
radial, and 18 times by 3 g radial in the two subjects 
tested (Table 2). The test was not attempted at 4 g, 
the definite impression being obtained that the barrier 
would have been insurmountable for these subjects at 
this acceleration. 

The results of these two tests show the tremendous 
increase in time required to accomplish simple locomo- 
tive movements of the body under the added stress of 
radial g. While the mat at the end of the centrifuge 
may have contributed to. the slowed movement of the 
subjects by offering considerable friction, nevertheless 
the inner surface of many combat aircraft would offer 
even greater frictional and obstructive resistance. 
When the subject was driven against the gymnasium 
mat by a force of 3 g, the time required to move the 
body 71/2 ft. (2.3 m.) or round a 22-in. (56-cm.) barrier 
was increased nine or more times. With a force of 4 g, 
most subjects were unable to,move their bodies. 

(3) Movement Against Radial g.—When testing 
movement against radial g, the subjects sat or stood 
against the end of the centrifuge and worked their way 
toward the center. Asa rule they were asked to walk, 
crawl, or pull themselves halfway to the center, a dis- 
tance of a little more than 81/2 ft. (2.6 m.). The nearer 
the subjects approached the center of the centrifuge, 
the less the radial g to which they were exposed. At the 
halfway mark the radial g was half that at the periphery. 
Under certain conditions similar circumstances might 
pertain in a spinning aircraft. The radial g on the 
centrifuge could have been kept constant by increasing 
the r.p.m. as the subject neared the center. No such 


TABLE 1 
Per Cent Increase in Time Required to Cross End of Centrifuge 
‘Under Radial g 


Time in 

Seconds 
Sub- atlg —————Per Cent Increase Under Radial g of ——— 
ject Gravity 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


2 1.16 148 272 383 601 690 cap 1,4 

12 177 276 397 «6518 843 1,735 Could not* 
211 1.62 920 Could not* 

Average 1.51 947 Three out 
of five 
could not 


* Subject could not get across end of centrifuge. 


TABLE 2 
Per Cent Increase in Time Required to Round Barrier 


Per Cent Increase Under 


Time in Seconds ial g of 
Subject at 1 g Gravity 1.0 2.0 3.0 
6 1.85 222 608 2,119 
8 1.79 240 620 1,453 
Average 1.82 231 614 1,791 


TABLE 3 


Effect of Radial g on Ability of Two Subjects to Walk Toward 
Center of Centrifuge 


Walked Distance of 81/2 Ft. 
With rope to 


Subject Radial g Unaided pull on 

2 0.5 Yes 
1.0 Yes aoe 
1.5 No Yes 
2.0 No Ves 
3.0 No | No 

8 0.5 Yes xs 
1.0 No Yes 
2.0 No No 

TABLE 4 


Effect of Radial g on Ability of Two Subjects to Crawl Toward 
Center of Centrifuge Unaided and on a Ladder 


Crawled Distance of 8'/2 Ft. 


Subject Radial g Unaided On ladder 

2 0.5 Yes 4 
1.0 Ves Yes 
1.5 No Yes 
2.0 No Yes 
2.5 No* 
3.0 Not 

8 0.5 Yes 
3. No 

* Was able to stretch out on ladder with feet on floor but could 
not get up even first rung. 


t Could not even get started on ladder. 


correction was made, however, in this study. The 
subjects wore rubber-soled shoes or overshoes which 
gave them some traction on the smooth plywood floor 
of the centrifuge. 

(a) In the first test the subjects simply walked or 
attempted to walk from the end of the centrifuge half- 
way to the center. When an acceleration had been 
reached at which they could not make any progress, 
they were given a rope fastened at the center of the 
centrifuge with which they could pull themselves to- 
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ROUNDING BARRIER DURING CENTRIPETAL ACCELERATION ' 
START SECONDS. | 10 SECONDS 


G20 


4.0 


1,465 
ould not* 


SECONDS 


91 23 
Toward 
Ft. 
ype to 
on 
oward Rounding a barrier during exposure to 2 g and 3 g radial. Notice time taken at 3 g. During control tests at 1 g gravity, 
only 1.8 sec. were required for the same maneuver (Table 2, subject 8). 
Ft. 
der 
ward the center. The test was designed to give some (b) In the second test the subjects were asked to 
measure of the difficulties confronting a man attempting crawl the 81/2 ft. (2.6 m.) toward the center, first with- 
to make his way to an escape hatch against radial g out assistance and then with the aid of a stepladder 
: and of how much a rope to pull on would help him. placed on the floor of the centrifuge. While doing this 
— Two subjects were tested. Up to 1 g radial their ability 
to cover the 8'/, ft. (2.6 m.) depended mainly on the 
ould traction between their shoes and the floor of the centri- 
fuge. Rubber soles were essential. The 1 g gravity 
aided them since it held them to the floor. Had gravity 
not been acting, the task would have been identical 
The with asking a man to walk up a vertical wall. At ac- 
tich celerations greater than 1 g the unaided subjects made 
oor no real progress; giving them a rope to pull helped only 
slightly (Table 3). Even with the rope only one of the 
or subjects was able to pull himself to the center at 2 g, 
alf- and at 3 g with the rope he failed. 
Pen At 2 g the test is not unlike requiring a man to pull 
»ss, himself 81/2 ft. (2.6 m.) up a rope with a man on his 
the back whose body weight is equal to his own, a difficult Fic. 4. At 4 g the subject is able to roll himself up the back of 


undertaking for most men. the centrifuge. 
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7'/2 ft. (2.3 m.). Five subjects were tested. On the 

average the time required for this maneuver was in- 
creased two and one-fourth times by 1 g radial, five times 
by 2 g radial, and nine times by 3 g radial (Table 1). 
Three of the five subjects tested could not accomplish 
the task at 4 g (Fig. 2). At or slightly above 4 g all 
of the subjects were immobilized. They could move 
their arms and legs, and these they pushed back and 
forth against the end of the centrifuge in an apparently 
aimless fashion. Their bodies were pinned to the end 
of the centrifuge with such force that they simply could 
not make any progress. While the gymnasium mat at 
the end of the centrifuge provided good “footing,” it 
also offered considerable friction to body movement. 
With suitable hand and foot grips some progress at 4 g 
might have been possible had a smoother, harder sur- 
face been provided. At 4 g, while forward or backward 
crawling was impossible, some of the subjects were 
able to roll themselves up the back of the centrifuge 
(Fig. 4). This motion might be of value in reaching an 
escape hatch. 

(b) Rounding a barrier. The motion required to get 
past the barrier was not unlike that which would be 
needed to get by a bulkhead in an aircraft (Fig. 3). 
The time required to round the barrier was increased 
two and one-fourth times by 1 g radial, six times by 2 g 
radial, and 18 times by 3 g radial in the two subjects 
tested (Table 2). The test was not attempted at 4 g, 
the definite impression being obtained that the barrier 
would have been insurmountable for these subjects at 
this acceleration. 

The results of these two tests show the tremendous 
increase in time required to accomplish simple locomo- 
tive movements of the body under the added stress of 
radial g. While the mat at the end of the centrifuge 
may have contributed to. the slowed movement of the 
subjects by offering considerable friction, nevertheless 
the inner surface of many combat aircraft would offer 
even greater frictional and obstructive resistance. 
When the subject was driven against the gymnasium 
mat by a force of 3 g, the time required to move the 
body 71/2 ft. (2.3 m.) or round a 22-in. (56-cm.) barrier 
was increased nine or more times. With a force of 4 g, 
most subjects were unable to,move their bodies. 

(3) Movement Against Radial g.—When testing 
movement against radial g, the subjects sat or stood 
against the end of the centrifuge and worked their way 
toward the center. Asa rule they were asked to walk, 
crawl, or pull themselves halfway to the center, a dis- 
tance of a little more than 81/2 ft. (2.6 m.). The nearer 
the subjects approached the center of the centrifuge, 
the less the radial g to which they were exposed. At the 
halfway mark the radial g was half that at the periphery. 
Under certain conditions similar circumstances might 
pertain in a spinning aircraft. The radial g on the 
centrifuge could have been kept constant by increasing 
the r.p.m. as the subject neared the center. No such 


TABLE 1 


Per Cent Increase in Time Required to Cross End of Centrifuge 
‘Under Radial g 


nds 

Sub- atlg Per Cent Increase Under Radial g ofp ———— 
ject Gravity 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


2 1.16 148 272 383 601 690 . 1,465 
6 1.85 1,081 Could not* 
1.13 177 276 397 «518 843 =:11,735 id not* 
211 1.62 920 Could not* 
Average 1.51 owes 947 Three out 


* Subject could not get across end of centrifuge. 


TABLE 2 
Per Cent Increase in Time Required to Round Barrier 


Per Cent Increase Under 


Time in Seconds Radial g of. 
Subject at 1 g Gravity 1.0 2.0 3.0 
6 1.85 222 608 2,119 
8 1.79 240 620 1,453 
Average 1.82 231 614 1,791 


TABLE 3 
Effect of Radial g on Ability of Two Subjects to Walk Toward 
Center of Centrifuge 


Walked Distance of 8'/2 Ft. 


With rope to 
Subject Radial g Unaided pull on 
2 0.5 Yes 
1.0 Yes 
1.5 No Yes 
2.0 No Yes 
3.0 No No 
8 0.5 Yes aay 
1.0 No Yes 
2.0 No No 
TABLE 4 


Effect of Radial g on Ability of Two Subjects to Crawl Toward 
Center of Centrifuge Unaided and on a Ladder 


Crawled Distance of 8'/2 Ft. 


Subject Radial g Unaided On ladder 

2 0.5 Yes Bes 
1.0 Yes Yes 
1.5 No Yes 
2.0 No Yes 
2.5 No* 
3.0 Not 

8 0.5 Yes 
1.0 No 

* Was able to stretch out on ladder with feet on floor but could 
not get up even first rung. 


t Could not even get started on ladder. 


correction was made, however, in this study. The 
subjects wore rubber-soled shoes or overshoes which 
gave them some traction on the smooth plywood floor 
of the centrifuge. 

(a) In the first test the subjects simply walked or 
attempted to walk from the end of the centrifuge half- 
way to the center. When an acceleration had been 
reached at which they could not make any progress, 
they were given a rope fastened at the center of the 
centrifuge with which they could pull themselves to- 
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LIMITING EFFECT OF CENTRIPETAL ACCELERATION 


ROUNDING BARRIER DURING CENTRIPETAL AGCELERATION | 
SECONDS: 10 ‘SECONDS 


RADIAL 
G 2.0 


ntrifuge 


START 


RADIAL 


uld not* 


ree out 
f five 
ould not 
1ake it 


ward the center. The test was designed to give some 
measure of the difficulties confronting a man attempting 
to make his way to an escape hatch against radial g 
and of how much a rope to pull on would help him. 
Two subjects were tested. Up to 1 g radial their ability 
to cover the 8/2 ft. (2.6 m.) depended mainly on the 
traction between their shoes and the floor of the centri- 
fuge. Rubber soles were essential. The 1 g gravity 
aided them since it held them to the floor. Had gravity 
not been acting, the task would have been identical 
with asking a man to walk up a vertical wall. At ac- 
celerations greater than 1 g the unaided subjects made 
no real progress; giving them a rope to pull helped only 
slightly (Table 3). Even with the rope only one of the 
subjects was able to pull himself to the center at 2 g, 
and at 3 g with the rope he failed. 

At 2 g the test is not unlike requiring a man to pull 
himself 81/2 ft. (2.6 m.) up a rope with a man on his 
back whose body weight is equal to his own, a difficult 
undertaking for most men. 


23 SECONDS 


Fic. 3. Rounding a barrier during exposure to 2 g and 3 g radial. Notice time taken at 3 g. During control tests at 1 g gravity, 
only 1.8 sec. were required for the same maneuver (Table 2, subject 8). 


(b) In the second test the subjects were asked to 
crawl the 8'/, ft. (2.6 m.) toward the center, first with- 
out assistance and then with the aid of a stepladder 
placed on the floor of the centrifuge. While doing this 


Fic. 4. At 4 g the subject is able to roll himself up the back of 
the centrifuge. 
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Fic. 5. Difficulty in donning a parachute during exposure to radial g. Control observations were made at 1 
time in seconds; notice added time necessary at 2 g and 3 g radial. 
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LIMITING EFFECT OF CENTRIPETAL ACCELERATION 


test the subjects wore electrically heated suits and back 
parachutes. Their performance while crawling was no 
better than that while walking (Table 4). Unaided at 
accelerations greater than | g they could not reach their 
objective. They reached it with the ladder at 2 g, but 
at 3 g they made no progress even with the ladder. 

(c) The ability of one subject (subject 2) to raise 
himself from a bomber seat under radial g was tested. 
The subject was able to get out of the seat and stand up 
at 2g; at 2.5 g with or without ropes to pull on he was 
able to raise himself in the seat but unable to get out 
of it and stand. At 3 g he could not raise himself from 
the seat. 

The experiments indicate what little progress man 
can make againgt radial g of any magnitude. In the two 
subjects tested, effective movement of the body while 
crawling, walking, climbing up a rope or ladder, or 
rising from a seat ceased between 2 to 3 g. Three g 
occurs at a radius of 20 ft. (6 m.) when the r.p.m. is 
only 21 and at a radius of 80 ft. (24 m.) when the r.p.m. 
is only 10.5 (Table 5). 

(4) Donning a Parachute When Exposed to Radial g. 
—The time required by three subjects to don.a standard 
back parachute was detérmined at 1 g gravity and at 
various increments of radial g. The average time re- 
quired at 1 g gravity was 17 sec. This was increased to 
21 sec. at 1 g radial, 41 sec. at 2 g radial, and 1 min. and 
15 sec. at 3 g radial—one of. the subjects failing to be 
able to don the parachute at all at this acceleration. It 
was obvious to the subjects and those observing the 
tests that at a little more than 3 g and certainly at 4 g 
none of them would have been able to get into the 
parachute (Fig. 5). 


COMMENT 


The series of subjects used in this study is small. 
The findings, however, are sufficiently striking and 
definite to show the extremely low accelerations at which 
aman’s efforts to move his body may become futile and 
to call attention to the serious limitations that increased 
gfrom any cause may place on a flier’s ability to escape 
from his aircraft. At an acceleration greater than 
2.5 or 3.0 g, or possibly even less, the flier may need 


assistance or conveniences not at present provided if he 


is to escape from his plane. 
It is possible that many aviators have a false sense 
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TABLE 5 
The R.P.M. at Various Radii Necessary to Give 3 g 
Radius in Feet R.P.M. 
10 29.7 
20 21.0 
40 14.8 
60 12.1 
80 10.5 
100 9.4 


of security under high acceleration because experience 
has taught them that in the ordinary upright sitting 
position they can move their arms and legs in the 
transverse direction at 6, 7, 8, or 9 g. Had they but 
attempted to raise their bodies or even their arms or 
legs against the force, they would have found how 
impossible it would have been for them to move from 
their seats. 

The study clearly shows that fliers will often need 
help if they are to escape from airplanes while increased 
accelerative forces are acting. The study emphasizes 
the need, recognized for some time by others, of addi- 
tional facilities and conveniences for escape from air- 
craft, and it indicates the importance of further con- 
sideration and development of such items as: (1) con- 
venient escape exits placed as pairs on opposite sides 
of the plane at all pilot and battle stations, (2) me- 
chanical or explosive opening of escape exits, and (3) 
ejection from the plane of the flier in a seat from his 
pilot, passenger, or battle station. 


CONCLUSIONS 


The study demonstrates that if fliers are to escape 
from aircraft when accelerations of 3 g or more have 
developed they will need help. It thus emphasizes the 
need, recognized for some time by others, of further 
consideration and development of devices designed to 
assist aviators when escaping from aircraft. 
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Jet Airplane Range Considerations 


JULIUS JONAS* 
Northrop Aircraft, Inc. 


ABSTRACT 


The lift coefficient for maximum range is both graphically and 
analytically determined. Engine overall efficiencies are discussed 
in their relation to range considerations. A maximum range ex- 
pression is derived, first for the case of constant altitude and then 
for the case of constant engine and air speed. The two expres- 
sions are compared as to their relative values. Climbing range 
and, finally, maximum overall range are discussed. The article 
is concluded by illustrating the range increase with altitude for a 
given amount of fuel, both by neglecting and by taking into 
account compressibility. 


SYMBOLS 


overall efficiency of conversion of fuel energy into 
thrust power 

= propulsive efficiency 

cycle efficiency 

airplane gross weight (also, initial gross weight), Ibs. 

available fuel load, Ibs. 

Q fuel consumption, lbs. per hour 

Y thrust specific fuel consumption, lbs. per hour per Ib. 

F net thrust, Ibs. 

S = wing area, sq.ft. 

R 

V 

V, 


= 


= 


= range, miles 
= freestream air speed, ft. per sec. 
gz = jet velocity, ft. per sec. 

Cy = fuel heat value, B.t.u. per Ib. 

a = mechanical equivalent of heat, 778.26 ft.lbs. per B.t.u. 
' M_ = jet mass air flow, slugs per sec. 

th = time to climb, sec. 

R.C. = rate of climb, ft. per sec. 

a, b, m, Ki, Kz, and C = constants defined in text 
Subscripts 

O = sea level 

h = altitude 

a = available 

R = required 

C = climb 

I = initial 

F = final 

L = cruise 

G = glide 

Symbols not otherwise identified are standard N.A.C.A. sym- 
bols. 


Lirt COEFFICIENT FOR MAXIMUM RANGE 


ee THRUST HORSEPOWER required for a jet airplane 
is given as in the case of the conventional airplane 
by 

t.hp.r = CpSpV*/1,100 (1) 


The thrust horsepower available may be written as 
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t.hp., = JQCyn,/550 (2) 
Combining Egs. (1) and (2) 
Cy = 2J Sp (3) 


Considering the cruising range of any airplane polar, 
the drag coefficient Cp and the lift coefficient C, may 
be conveniently expressed as a straight line variation, 
(See Fig. 1.) Thus C, may be written as 


CG, = aCp —b 
or 
Cy = (G,/a) + (6/a) 


Remembering the basic lift equation, Cp may be 
written as 


Cp = (2W/aSpV*) + (b/a) (4) 
Substituting Eq. (4) into Eq. (3) and solving for Q 
Q = (WV/aJCyn,) + 
Indicating by 


Ki, = W/aJCy (6) 
and by 
Ky = bpS/2JCya (7) 
Eq. (5) may be written as | 
Q = (1/m,)(KiV + KeV*) (5’) 


Dividing both sides of Eq. (5’) by V and differentiating 


CL 


AIRPLANE POLAR 


SLOPE a 


Linear representation of airplane polar cruising range. 
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dQ/V)_ dm 1 
+ K2V?) av + 


At the speed for which Q/ V becomes a minimum, 
d(Q/V)/dV = 0 
Rearranging 
dy,/dV = ,[2K2V/(Ki + K2V)] (8) 
Integrating 
dni/n = S [2K2V/(Ki + K2V2)]dV 


Solving 
nm = C(Ki + K2V?) (9) 


Varying the value of the integration constant C in 
Eq. (9), a family of curves is obtained. The constant C 
is defined by the boundary condition that Eq. (9) must 
have a common point of tangency with the actual rela- 
tionship of overall efficiency vs. air speed. 


OVERALL EFFICIENCY 


' The overall efficiency of a jet engine is given by the 
ratio of thrust horsepower developed to the heat energy 
input. 


The above expression can also be written as the 
product of the cycle efficiency 7, and the propulsive 
efficiency Thus, 


= Mp 


Whereas the thermal efficiency is a measure of how 
much of the original fuel energy is being converted into 
kinetic energy, the propulsive efficiency indicates how 
much of the kinetic energy becomes thrust power. 

The overall efficiency will normally tend to remain 
essentially constant over an appreciable range of 
engine operating speeds, provided ambient temperature 
and air speed is maintained constant. 

The range of engine speeds operating at constant 
overall efficiency corresponds to the peak of the overall 
efficiency curve and thus represents the best economy 
range of engine speeds. Beyond this range operating 
speeds will reach maximum thrust values with, for the 
most part, somewhat decreased overall efficiencies. 
The maximum allowable speed giving maximum thrust 
may thus be termed best power speed. 

Obviously engine speeds below the range of best 
economy overall efficiencies are of little interest. The 


overall efficiency will be, usually, a uniquely determined 
straight line function of the air speed within the above- 
mentioned range of best economy engine speeds. 

The overall efficiency will increase in most cases with 
altitude at a given engine speed and true air speed. 
This is brought about by a decrease in ambient air 
temperature. 


In fact, variations in overall efficiency 


become negligible above the lower level of the iso- 
thermal atmosphere. 


MAXIMUM RANGE ANGLE OF ATTACK 


In most instances it will be possible to express the air 
speed-overall efficiency relationship in a linear form by 


where m is of the order of magnitude of approximately 
5,000 to 6,000. The boundary conditions will thus be 
satisfied by putting in Eq. (8) 

dn,/dV = 1/m 


and 
n = V/m 


Solving for V 
V=V Ki/K2 (11) 


Substituting the original values for K, and K2 as defined 
in Eqs. (6) and (7), Eq. (11) becomes 


V = V2W/bpS (11’) 


Thus the lift coefficient to be maintained for maxi- 
mum range becomes simply equal to d. 

This suggests a simple method for the determination 
of b. (See Fig. 2.) It consists of placing a straight 
edge tangent to the polar and changing the point of 
tangency until C, becomes equal to b. By the same 
token the value of C,/Cp for maximum range becomes. 
equal to a/2. 


CL 
AIRPLANE ATTITUDE 
FOR MAXIMUM RANGE 
b 
b 


Fic. 2. Determination of lift coefficient for maximum range. 
Linear relationship of air speed vs. overall efficiency. 


Considering the general form of an airplane polar 
Cp = f(G) 


The above-determined condition for maximum range 
may also be expressed as 


dCp/dC,, = 2C,/Cp 


V = mn (10) 
(2) 
(3) | 
, May 
| 
7 
) 
(5) 
(6) 
(5’) 
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Rearranging and dividing by C,'/* 
— C,*(dCp/dC,) = 0 
This may also be written as 
= 0 
and clearly represents the condition of 


= maximum 


Thus the above method is merely a graphical solution 


for the maximum value of C,'” yg 


MaxIMuUM RANGE AT A GIVEN CONSTANT ALTITUDE 


If constant-altitude maximum-range cruising is re- 
quired, it will be necessary to maintain a constant angle 
of attack giving the lift-drag ratio for maximum range 
as determined above. This is accomplished by de- 
creasing engine speed and flight speed as the fuel load is 
consumed. It is assumed that the engine speed range 
covered lies on the constant value portions of the overall 
efficiency curve. 

The incremental range dR may be written as 


dR = Vdt 
Also 
dt = (V/Q)dW 
Thus 
dR = (V/Q)dW 


V may be written as 


V = WV 2/C,Sp 


and Q 
Q = 
Substituting 
it 2/SCp W-"*dW 
m Cp 


Integrating between W and W — Wy and using consis- 
tent units 
Cy 


5,280m Co 2/SCip [W (W — 


(12) 


Substituting for C,/Cp = a/2 and introducing the initial 


and final speed V; and V;, respectively, 
R = (aJCy/5,280m)(Vr — Vr) (13) 


Sometimes it may be more convenient to use the thrust 
specific fuel consumption + (Ibs. per hour per Ib.) 


y = 3,600m/JCy 
Thus 
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Maximum RANGE AT CONSTANT ENGINE SPEED AND 
AIR SPEED 


If an airplane is originating its range flight at a given 
altitude with a given initial gross weight, the angle of 
attack for maximum range may again be selected, 
This will define both the engine speed and air speed re- 
quired. If the initial angle of attack and engine speed 
are maintained constant, the airplane will slowly climb 
maintaining constant true air speed. 


At the initial point 


2W/p,SC, 


At the final point 
Vp = V2(W — Wr)/prSCr 


Maintaining constant air speed and angle of attack, the 
following relation is given between gross weight and 
altitude density 


(W = Wr)/W 


However, in maintaining constant engine speed and 
angle of attack, the air speed will remain constant only 
if the initial and final thrust values vary as o. 


Fy / Fy 


(15) 


(16) 


This condition is met, however, since any jet airplane 
range flight is invariably calculated to occur above the 
lower level of the isothermal atmosphere. 

It can be shown that for constant ambient tempera- 
ture, constant engine and air speed, both propulsive and 
cycle efficiency will remain constant. Thus, both 
engine weight flow and thrust will vary as co. 


dR = (V/Q)dW 
dW/de = 


F,V/Q; = FV/Q 
Q 


Substituting 


dR = W(V/Q:)(do/c) 
R = W(V/Q,) log. (¢1/¢r) 


Substituting for 


~ 1.4670, W— We 


Thus, for the initial selected altitude V may be calcu- 
lated for the gross weight W and the maximum range 
lift coefficient b. Hence, V/Q may be obtained from 
available engine data. 

Alternate expressions may be written in terms of the 
overall efficiency slope m or the thrust specific fuel 
consumption 


(17) 


aVJCy W 
W, 


= — 
10,560m W F 
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(18) 


2.934 
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COMPARISON OF MaximuM RANGE FOR CONSTANT 
SPEED WITH THAT FOR CONSTANT ALTITUDE 


Taking the ratio of the maximum range when flying 
at constant airplane and engine speed as well as at 
maximum range lift coefficient—R,—to the maximum 
range when flying at constant altitude and at the max- 
imum range lift coefficient—Rz,—1t easily follows from 
Eqs. (14) and (19) 


_ 
Rin 


log, 1/[1 — (Wr/W)] 
21 — VY 1— W,/W) 


For the trivial case of W/W approaching zero, the 
ratio R,/R,, approaches unity. At a representative 
value of Wy/W = 0.3, R, is approximately 8 per cent 
larger than R,,. Actually, thé ratio R,/Rz, will 
always be larger than indicated by Eq. (20). The in- 
crement is due to the potential gliding distance ARg 
from the final altitude to the initial altitude. 


ARg = Cp) max. Ah 


where Ah is defined by o; and ar. 

Again taking a representative case of W/W = 0.3 
and (C,/Cp) max. = 12. Ah becomes equal to approx- 
imately 7,500 ft. 

Thus, AR, = 17 miles, which is only a small fraction 
of the total range. 


(20) 


RANGE IN CLIMB 


The average value of V;/Q for a climb to a given 


altitude h is given by 
ve) - (© 
== dh 
( 


By the same token the actual horizontal range covered 
in the course of the climb is given by 


th (RS). 
22 


Thus, the fuel required to climb to altitude, Weg, is 


given by 
=) ah 
1467 V, R. 
Le 
Ve 


It is obvious that in order to obtain maximum overall 
range, the fuel expended in the course of the climb to a 
given altitude must become a minimum as the fuel is 
more efficiently utilized in altitude maximum-range 
level flight. 


(21) 


Ren 


Wre = (23) 
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If the available thrust of a given airplane were to be 
increased by a given amount, as, for instance, by in- 
creasing the engine speed at a given airplane speed, the 
following effects would take place at any given alti- 
tude: 

(a) The overall efficiency would remain essentially 
the same. 

(b) The V;-/Q value would vary inversely with the 
thrust and thus would be subject to essentially the same 
percentage change in absolute value. 

(c) The excess thrust horsepower would increase by 
a larger percentage than the thrust horsepower avail- 
able. The time required to climb also would corre- 
spondingly decrease by a larger percentage. 

Thus it appears that minimum climb fuel consump- 
tions, Wr,, and correspondingly the maximum overall 
range will be obtained by climbing with maximum 
available thrust. This reasoning is based on the 
assumptions that both overall efficiencies and best climb 
speeds remain essentially unaffected by a maximum 
thrust increase. However, as a rule it will still hold 
true for appreciable deviations from these assumptions 
and within which most practical cases will fall. In any 
event, the above expression for Wr, must become a 
minimum. 


MAXIMUM OVERALL RANGE 


The range is composed of a climb range Re to altitude 
h, a constant engine and air-speed range R, with an in- 
herent altitude increase until all fuel is exhausted, and a 
gliding range R, from altitude h + Ah. 

Thus 


(24) 


R= Re + + Rg 
The altitude / selected which will give maximum overall 
range is given by 


dR,/dh = (25) 


—d(Re + Rg)/dh 
This may be interpreted as follows: 

With increasing altitudes less and less fuel becomes 
available for the range R, that may be attained. 

Miles-per-gallon values increase with increasing 
altitudes. 

Clearly, there will be an altitude at which R,, will 
become a maximum and beyond which R, will decrease. 
On the other hand (Re + Rg) will increase with altitude 
until all fuel is exhausted. 

Thus 


d(Re + Rg)/dh 


will always be positive. dR,/dh is positive up to the 
altitude of Rymez.. Above that altitude it becomes 
negative. 

It is interesting to note, therefore, that maximum 
overall range will result at an altitude somewhat higher 
than that at which R,; becomes a maximum. 


— Wr 
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RANGE INCREASE WITH ALTITUDE FOR A GIVEN AMOUNT 
OF FUEL 


In studying the effect of altitude on range, it is 
assumed that a linear relationship exists between overall 
efficiency and air speed. As a first approximation, let 
it be assumed that the compressibility effect may be 
neglected. Thus, the sea-level airplane polar may be 
used for all altitudes. In assuming a unique polar, the 
same maximum-range angle of attack will apply for all 


altitudes. 

Thus, 

Va = (26) 
and 
t.hp.tg = 
Further 
t.hp.z, = 

Thus 


Substituting for », in Eq. (27) and remembering Eq. 
(26), Eq. (27) becomes 


(Q/V)n/(Q/V)o = (28) 


Therefore, at approximately 40,000 ft. altitude the jet 
airplane range would be twice that at sea level, assuming 
m,/m) = 1.0. An actual representative value of m,/mo 
is about 0.9, thus giving a still larger increase in range 
with altitude. 


Macu NuMBER EFFECT 


Under actual conditions, increasing the altitude will 
result in a modification of the airplane polar. Gen- 
erally, the minimum drag coefficient will occur at a 
higher lift coefficient than at sea level. Normally, high- 
speed airplanes, such as fighter and pursuit airplanes, 
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will have a sufficiently high critical Mach Number 
keep the value for maximum range out of the comm 
pressibility range. In the case of the bomber and carg@ 
airplane, this will not normally be true. In the latte 
case it may be found that the angle of attack at altitud 
required for maximum range will increase compared t 
sea level, resulting in a larger lift coefficient. Curiously 
enough, the lift-drag ratio tends to remain substantially 


constant. On that basis 


Va = 
and 
t.hp.ig = 
Finally, Eq. (29) may be written similarly to Bal 
(28) as 
vy, (m\f 


Thus the range at altitude is approximately reduced by§ 
the ratio of Czo/Cz, as compared to the case where the 
maximum range lift coefficient falls outside the coma 
pressibility region at all altitudes. 

This emphasizes the importance of striving for high 
critical Mach Numbers, not only from the standpoint} 


of attaining good high-speed characteristics, but also taj 


insure optimum range performance. More knowledge 


is required also about functional jet engine operationj™ 


since this seems to be a conditio sine qua no for the | 
attainment of high overall efficiencies. 
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